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INTRODUCTION 


This  section  includes  a  discussion  of  the  objectives  of  the  research  pro¬ 
gram,  a  review  of  the  work  performed,  and  the  results  obtained.  The  work 
was  performed  by  The  Ohio  State  University  from  15  July  1973  to  30  June 
1974. 


OBJECTIVE 

The  objective  of  this  program  was  to  investigate  the  dynamic  effects  en¬ 
countered  at  the  gear  mesh  due  to  incipient  gear  tooth  failure  modes. 

The  specific  failure  mode  under  study  was  pitch-line  pitting  of  a  gear 
tooth. 

A  goal  of  the  study  was  to  develop  a  dynamic  model  capable  of  predicting 
gear  behavior.  Experimental  data  were  to  be  used  in  verification  of  the 
model  and  the  results  then  applied  to  possible  applications  in  the  evalua¬ 
tion  and  development  of  gear  diagnostic  techniques. 


PROGRAM  CONDUCT 


The  program  was  divided  into  five  tasks: 

Task  I.  Mathematical  Formulation 

Task  II.  Model  Simulation 

Task  III.  Experimental  Studies 

Task  IV.  Model  Verification 

Task  V.  Recommendations 


During  Task  I  a  set  of  equations,  based  on  analytical  considerations,  was 
developed  to  predict  the  dynamic  response  of  a  geared  system  both  when  the 
system  contains  all  good  gears  and  when  gear  tooth  faults  are  present.  The 
basic  procedure  was  an  analysis  of  the  torsional  dynamic  response  of  the 
system  with  and  without  gear  tooth  faults.  A  separate  formulation  was  per¬ 
formed  to  determine  gear  tooth  stiffness  functions. 

The  work  of  Task  II  consisted  of  using  the  equations  developed  during  Task 
I  to  develop  a  computational  procedure  for  evaluating  the  geared  system’s 
response  characteristics.  The  simulation  was  performed  on  The  Ohio  State 
IBM  System  370/165  digital  computer.  Programming  was  performed  in  FORTRAN 
and  utilized  a  number  of  IBM  subroutines  and  subprograms. 

A  test  apparatus  was  designed  and  built  and  experiments  were  performed  in 
Task  III.  The  test  stand  utilized  a  Generator  Drive  Offset  Quill  gearbox 
which  is  used  on  the  UH-1  helicopter.  Faults  were  manufactured  into  a 
total  of  four  gears;  experiments  were  performed  on  good  gears  as  well  as 
those  containing  the  faults.  A  number  of  parameters  were  monitored  during 
testing,  including  acceleration  at  various  gearbox  positions  and  dynamic 
strain  of  the  gear  mesh. 
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The  results  of  the  simulation  and  experiments  were  compared  as  part  of 
Task  IV*  Appropriate  adjustments  in  the  model  were  made  as  a  result  of 
experimental  evidence* 

Task  V  consisted  of  evaluating  the  results  of  the  simulation  with  respect 
to  possible  applications*  These  consist  of  the  applicability  of  the  model 
and  modelling  procedures  to  the  evaluation  and  development  of  diagnostic 
procedures ,  as  well  as  further  developmental  possibilities  of  the  simula¬ 
tion  which  might  prove  valuable. 


REPORT  CONTENTS 


The  testing  equipment  is  described  and  a  theoretical  analysis  performed  on 
the  test  system.  Data  are  presented  on  time-domain  response  and  frequency- 
domain  response  from  both  experiments  and  simulation  results.  The  data  are 
presented  for  both  the  cases  of  operation  of  the  system  in  good  condition 
and  also  when  faults  are  present.  Changes  noted  as  the  result  of  fault 
introduction  are  described,  and  correlation  between  theory  and  experiment 
is  discussed. 
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BACKGROUND 


This  study  was  performed  to  develop  a  dynamic  model  of  a  geared  system 
which  incorporates  the  effects  of  gear  tooth  surface  faults.  The  response 
characteristics  of  the  model  could  then  be  used  to  determine  differences 
between  the  vibrations  generated  by  good  gears  and  those  possessing  faults. 
Once  such  a  capability  is  available,  it  may  be  applied  to  the  development 
and  evaluation  of  discriminants  for  use  in  vibration  diagnostic  procedures. 

A  model  of  this  sort  must,  of  course,  be  verified  as  reliable  through  cor¬ 
relation  with  experimental  data.  However,  once  such  a  me  1  is  considered 
valid,  variations  in  the  system  such  as  load,  speed,  fault  .  *.e,  fault  posi¬ 
tion,  etc.,  may  be  studied  without  the  immediate  need  for  experimental  data. 
This  type  of  result  can  present  both  cost  and  time  savings  in  the  develop¬ 
ment  and  analysis  of  diagnostic  techniques. 


DIAGNOSTIC  ANALYSIS 


A  great  deal  of  interest  has  been  evidenced  in  the  use  of  vibration  diag¬ 
nostic  techniques  in  conjunction  with  maintenance  procedures.  Diagnostic 
techniques  are  defined  as  those  which  identify  the  condition  of  a  system 
from  its  external  signs  or  symptoms.  In  the  case  dealt  with  in  this  study, 
the  condition  of  a  power  train  gear  is  to  be  determined  by  analyzing  vibra¬ 
tion  signals  from  the  power  train  housing.  In  actual  practice  such  a  pro¬ 
cedure  eliminates  the  need  for  visual  inspection  of  the  individual  power 
train  component.  There  has  been  considerable  study  of  vibration  diagnostic 
procedures  as  applied  to  power  trains.  An  overview  may  be  obtained  from  a 
report  on  the  subject  prepared  by  Houser  and  Drosjack.1  A  number  of  tech¬ 
niques  have  shown  strong  promise  for  application  with  bearing  failures. 
However,  the  same  statement  may  not  be  made  with  respect  tc  gearing. 

In  bearing  diagnostics,  the  techniques  which  have  shown  the  most  success 
analyze  in  one  manner  or  another  the  impact  which  a  rolling-element  surface 
fault  produces  as  the  fault  passes  through  the  interface  between  two  con¬ 
tacting  surfaces.  This  force  then  excites  bearing  support  structures  which 
produce  a  measurable  "ringing"  response  as  seen  in  Figure  1.  The  repetition 
rate  relates  to  a  characteristic  frequency  which  may  be  determined  on  the 
basis  of  fault  position  and  is  readily  calculable.  The  "ringing"  frequency 
is  not  yet  predictable  theoretically. 

A  wide  variety  of  analysis  procedures  are  used  in  processing  the  vibration 
signals.  SKF,  in  the  Shock  Pulse  Meter,  uses  the  wave  propagation  response 
to  the  impact  force  to  "ring"  an  accelerometer  (mounted  on  the  machine 
structure)  at  its  resonance.2  The  output  of  this  accelerometer  is  monitored 
for  both  peak  amplitude  and  frequency  of  occurrence  of  rings.  Mechanical 
Technology  Incorporated,  with  their  Bearing  "Ring"  Analyzer  processes  the 
signal  as  shown  in  Figure  l,3  They  analyze  it  for  both  the  amplitude  of  the 
"ringing"  response  and  the  repetition  rate  of  the  impact.  General  Electric 
Company  uses  a  simpler  analysis  in  their  Impact  Index  Analyzer,4  where  they 
calculate  the  value  of  the  ratio  of  peak  vibration  signal  to  the  r.m.s. 
level.  These  are  just  three  examples  from  many  analysis  techniques.  .Almost 
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all  of  them,  however,  In  some  manner  make  use  of  information  relating  re¬ 
sponse  to  the  repetition  rate  of  the  impact  and  the  effect  of  fault  size, 
position,  and  type  on  amplitudes  of  bearing  force  and,  hence,  vibration 
that  may  be  expected. 

There  is  much  less  understanding  of  the  mechanisms  encountered  in  gearing 
when  a  fault  is  present.  The  dynamic  phenomena  have  not  been  analyzed  to 
any  great  extent,  and  data  obtained  experimentally  are  difficult  to  inter¬ 
pret.  A  fault  may  be  expected  to  cause  dynamic  forces  in  the  gear  mesh. 
However,  unless  sliprings  or  telemetry  devices  are  used,  this  force  cannot 
be  monitored.  Instead,  vibrations  on  housing  structures  of  the  gearbox  are 
usually  measured.  The  characteristics  of  the  transmission  path  between  the 
gear  mesh  and  the  transducer  may  severely  distort  the  signal  and  make  anal¬ 
ysis  of  it  quite  difficult. 

There  have  been  some  gear  diagnostic  attempts  based  on  the  observed  changes 
in  the  frequency  spectrum.  Faulty  gears  typically  exhibit  more  energy  at 
harmonics  of  the  meshing  frequency  and  at  sidebands  of  both  the  meshing 
frequency  and  its  harmonics.  The  interrelationships  between  position  and/ 
or  size  of  these  spectrum  changes  to  the  fault  type  and  size  have  not  been 
well  understood. 

General  Electric  Company  used  three  ratios  as  gear  condition  indicators  on 
their  TF-34  Analyzer. b  They  took  a  ratio  in  the  manner  of  their  Impact  Index, 
a  ratio  of  harmonic  content  to  the  amplitude  of  the  gear  meshing  frequency, 
and  a  ratio  of  side-band  content  to  the  fundamental  and  harmonic  content. 

A  joint  technique  was  developed  by  General  Electric  and  Mechanical  Technology 
Incorporated  for  monitoring  gear  wear.6*7  The  ratios  of  energy  at  the  first, 
second,  and  third  harmonics  to  the  fundamental  gear  meshing  frequency  were 
formed.  Experimental  results  showed  these  ratios  to  increase  as  the  gear 
wear  progressed.  Some  later  work  on  this  technique  included  preconditioning 
of  the  data  by  time-averaging  of  the  signal  at  the  gear-meshing  period  or 
shaft  rotational  period.8  Time-averaging  adds  together  blocks  of  data  at 
intervals  of  the  gear-meshing  period.  Response  which  is  uncorrelated  at 
this  period  will  cancel  itself  out,  thus  allowing  the  response  of  the  gear 
to  predominate.  No  widespread  application  of  any  of  these  gear  diagnostic 
techniques  has  been  discovered  by  the  investigators  in  this  study. 

One  of  the  primary  roadblocks  to  further  development  of  gear  diagnostic 
techniques  has  been  the  lack  of  analytical  results  relating  gear  response 
to  faults.  This  is  not  to  say  that  dynamic  studies  of  gearing  have  not 
been  preformed.  There  have  been  a  large  number  or  studies  performed,  but 
most  of  these  deal  with  either  design  procedures  of  the  dynamic  effects  of 
manufacturing  errors.  A  listing  of  work  would  include  study  of  dynamic 
loading,  mesh  stiffnesses,  load  distributions,  gear  noise,  and  geared  sys¬ 
tem  response. 


GEAR  DYNAMIC  LOADING 

An  area  where  interest  has  long  been  evident  is  the  study  of  acceptable 
load  rating  of  gear  sets.  In  the  nineteenth  century  a  number  of  theories 
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yielding  wide  variations  in  calculated  allowed  loads  were  used.  A  presenta¬ 
tion  was  made  in  1892  by  W.  Lewis  to  the  Engineers  Club  of  Philadelphia 
which  included  the  following  equation:8 


Fs  -  TdWPcY  (1) 


where  Fs  =  safe  allowable  load,  lb 
xd  «  design  stress,  lb /in. 2 
W  ■  facewidth  of  gear,  in. 

Y  «  form  factor  depending  upon  tooth  shape 
Pc  =  circular  pitch,  in. 


This  equation  was  based  on  an  analysis  of  the  material  working  stress  using 
simple  beam  theory.  Lewis  also  felt  that  gear  velocities  should  be  included 
in  the  equation.  This  led  to  the  development  of  the  following  equation  by 
Barth: 1U 


d  "  T  500 ~  (2) 

where  r  *  safe  static  stress,  lb/in.2 

v  =  pitch-line  velocity,  ft/min 

This  modified  design  stress  was  then  used  in  Lewis'  design  equations.  It 
was  one  of  the  earliest  dynamic  increment  or  derating  factors  used  to  ac¬ 
count  for  dynamic  loading  encountered  by  normal  gear  operation. 

In  1913>  E.  Buckingham  authored  a  report  of  the  A.S.M.E.  Research  Committee 
on  Dynamic  Loads11  which  presented  the  results  of  work  carried  on  by  the 
committee  on  the  effect  of  manufacturing  errors  and  pitch-line  velocity  on 
the  dynamic  loads  of  gears.  The  equations  which  were  developed  are  still 
in  use  with  appropriate  modifications. 

In  the  early  1950's,  Tuplin  presented  a  more  refined  analysis  for  deter¬ 
mining  the  loads  due  to  manufacturing  errors.12  Figure  2  shows  the  model 
used  in  his  analysis.  The  mass  (m)  represents  the  equivalent  gear  inertias, 
and  the  stiffness  (K)  is  an  average  mesh  stiffness  calculated  using  a 
static  analysis.  The  wedge,  moving  at  a  velocity  V,  represents  a  positive 
pitch  error  of  size  e.  Attia  published  the  result  of  experimental  work 
which  used  strain  gages  to  record  the  deflections  of  gear  teeth  and  the 
maximum  instantaneous  tooth  loading  as  the  gear  teeth  passed  through  the 
zone  of  contact  at  various  speeds  and  static  loads.13  His  results  showed 
differences  from  the  predictions  o°  both  Buckingham  and  Tuplin.  In  195*4, 
Reswick  presented  a  more  sophisticated  analysis  procedure  for  determining 
dynamic  loading.14  In  it,  he  used  a  parabolic  cam  to  represent  tooth 
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errors  and  allowed  the  inclusion  of  the  effect  of  more  than  one  tooth  pair. 
Figure  3  shows  the  model  which  was  used. 

IJlemann  and  Rettig  reported  the  results  of  a  study  which  experimentally 
determined  the  effect  of  various  operating  conditions  and  tooth  errors  on 
the  dynamic  loading  and  deflection  of  the  gear  teeth.15  Experimental  work 
was  performed  by  Harris  to  study  tooth  stresses  by  photoelastic  methods.13 
He  concluded  that  dynamic  loading  may  be  attributed  to  manufacturing  errors, 
parametric  excitation  due  to  stiffness  variations  in  the  mesh,  and  backlash 
effects  in  the  gears,  Gregory*  Harris,  and  Munro  performed  studies  which 
concluded  that  damping  was  an  important  consideration  in  any  analysis.17 
Later  work  by  Munro  stated  that  transients  did  not  die  out  as  quickly  as 
previously  considered,18  For  this  reason,  he  suggested  that  analyses  must 
consider  several  tooth  pairs  to  be  accurate,  Richardson  performed  a  com¬ 
prehensive  study  of  tooth  deflection,  stress,  and  loading  that  included 
both  theoretical  development  and  experimental on;19  his  model  is  shown 
in  Figure  4, 

A  planetary  geared  system  was  analyzed  by  Kasuba,20  who  developed  a  model 
as  shown  in  Figure  5.  The  errors  were  simulated  by  a  displacement  strip 
[S(t)  ] .  He  also  extended  his  model  to  two  degrees  of  freedom.  Hahn  later 
extended  Kasuba's  analysis  to  include  the  effect  of  other  system  components, 
including  shafting,  bearing,  and  housing.21 

Houser  performed  a  study  of  spur  and  helical  gears  which  led  to  the  devel¬ 
opment  of  a  semiempirical  dynamic  loading  factor.  »  3  He  calculated 
velocity  differences  based  on  pitch  errors  and  performed  experimental 
studies  to  determine  the  effects  of  pitch  errors  and  face-width  variations 
on  dynamic  tooth  loads.  Dynamic  strains  were  measured  with  the  use  of 
strain  gages  mounted  in  the  tooth  roots. 

Various  studies  of  a  nature  similar  to  many  previously  mentioned  have  also 
been  performed  by  a  number  of  Japanese  investigators.24”27 

In  the  analyses  discussed  in  the  above  paragraphs,  the  prime  considerations 
involved  either  normal  gear  operation  or  operation  with  gears  containing 
manufacturing  errors.  Most  of  the  work  also  analyzes  only  the  interaction 
between  one  tooth  pair  or  at  most  two  or  three  tooth  pairs.  While  the  in¬ 
formation  is  of  some  value  in  studying  the  effects  of  tooth  surface  faults 
on  gear  system  response,  the  studies  are  not  adequate  to  provide  the  de¬ 
sired  predictive  capabilities. 


MESH  STIFFNESS  AND  TOOTH  DEFLECTION  ANALYSIS 


A  number  of  studies  have  been  performed  on  gear  mesh  stiffness  functions 
and  deflections  under  load.  This  type  of  information  has  two  uses;  i.e., 
information  on  torsional  resonant  response  characteristics  may  be  obtained, 
and  the  values  of  deflection  may  be  used  to  modify  teeth  for  smoother 
operation  under  load. 
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Figure  4.  Richardson'*  Model  for  Gear  Analysis. 


Figure  5.  Kasuba’s  Single-Degree-of-Freedosn  Model. 


Work  in  this  area  considers  compliance  (the  inverse  of  stiffness)  as  well 
as  stiffness  itself.  The  mesh  compliance  coefficient  (C^)  is  defined  to 
be  the  deflection  in  inches  (measured  relative  to  the  base  circle  of  the 
gear)  per  pound  of  force  applied  along  the  pressure  line  per  inch  of  face- 
width  of  the  gear.  Typical  values  of  compliance  coefficient  versus  dis¬ 
tance  along  the  line  of  action  f  shown  in  Figure  6. 

Load  sharing  between  tooth  pairs  must  also  be  included  when  determining 
total  mesh  stiffness.  In  a  spur  gear  pair  having  a  contact  ratio  between 
one  and  two,  the  load  will  be  carried  first  by  one  pair,  then  by  two  pairs, 
then  by  one  pair,  etc.  This  change  in  loading  can  cause  abrupt  variations 
in  the  overall  mesh  stiffness.  Meshing  is  assumed  to  occur  only  along  the 
line  of  action  of  the  gears,  and  load-sharing  occurs  on  portions  of  this 
line  which  are  dependent  on  the  tooth  and  mesh  design  parameters,  and  which 
may  be  determined  geometrically.  When  more  than  one  tooth  pair  is  in  con¬ 
tact,  the  tooth  pairs  are  treated  as  springs  in  parallel.  A  second  possible 
analysis  procedure  which,  to  the  knowledge  of  the  investigators,  has  not 
been  performed  is  to  model  noninvolute  action  due  to  errors  and  tooth  de¬ 
flections.  This  would  move  contact  off  the  line  of  action  and  would  re¬ 
quire  an  iterative  computation  procedure. 

In  1926,  Timoshenko  and  Baud28  presented  an  equation  for  tooth  bending 
deflection  which  was  based  on  the  analysis  of  a  tooth  modelled  as  a  trap¬ 
ezoidal  beam  (see  Figure  7).  An  average  value  was  used  for  shear  def lec- 
tion^  and  the  equations  for  contact  between  two  cylinders  gave  Hertzian 
deflection  values.  The  equation  was: 


5=5+  5. 


12P  A  a  W  a 

=  E  h7  [\2  2l)  \i 


+ 


4p(  i-a)  (l-[j ) 
(h+h0)E 


+ 


(3) 


where  5  =  total  deflection 

5g  -  bending  deflection 
5S  =  shear  deflection 
5h  -  Hertzian  deflection 

rg,rp  =  radius  of  curvature  of  contacting  involute 
surfaces  at  contact 
b  -  width  of  contact  strip 
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TH  PAIR  COMPLIANCE,  in./lb  force /in. (xIO6) 


DISTANCE  FROM  PITCH  POINT  AUONG  LINE  OF  ACTION 
TIMES  DIAMETRAL  PITCH  (in.) 


Figure  6.  Typical  Plots  of  Tooth  Pair  Compliance  Versus 
Contact  Position* 
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Maximum  deflections  tinder  load  were  shown  to  be  negligible  in  comparison 
to  manufacturing  errors  of  the  day. 

In  a  subsequent  paper  publish?  I  in  1929 »  Baud  and  Peterson  used  Timoshenko's 
equations  to  obtain  the  variation  of  tooth  pair  compliance  as  a  function  of 
position  of  the  point  of  contact  along  the  arc  of  action.29  The  model  they 
used  for  tooth  bending  was  similar  to  Timoshenko's.  The  load  used  in  their 
equations  was  equal  to  the  component  of  tooth  force  which  is  perpendicular 
to  the  center  of  the  tooth,  with  the  deflections  being  geometrically  re¬ 
solved  to  the  pressure  line.  Figure  8  shows  the  total  compliance  and  the 
individual  compliance  components  as  calculated  for  a  single  tooth  pair. 
Figure  9  shows  the  mesh  compliance  when  load-sharing  was  included  in  their 
calculations . 

In  1938*  H.  Walker  performed  an  experimental  investigation  of  tooth  pair 
compliance.30*31  His  testing  apparatus  consisted  of  a  nonrotating  gear  to 
which  a  load  could  be  applied  by  a  lever  arm  arrangement.  Load  geometry 
could  be  changed  by  varying  the  position  of  the  gear  center  with  respect 
to  the  point  of  force  application.  Figure  10  shows  results  obtained  when 
the  load  at  a  given  contact  position  was  varied.  Figure  11  shows  the  re¬ 
sults  of  varying  contact  position  for  given  loads.  Based  on  his  experi¬ 
mental  data,  Walker  was  able  to  conclude  that  the  mesh  compliance  was  in¬ 
dependent  of  diametral  pitch.  He  also  showed  that  mesh  compliance  was 
essentially  a  linear  function  of  load.  Other  investigators  have  also  found 
these  two  ideas  to  be  true. 

A  more  rigorous  theoretical  approach  was  taken  by  Weber.32  His  analysis 
was  utilized  in  the  modelling  procedures  used  in  the  research  program  dis¬ 
cussed  in  this  report.  His  method  is  detailed  in  Chapter  4. 


The  above  procedures  determine  gear  tooth  deflections  on  the  basis  of  beam 
analysis.  In  the  case  of  wide  face-width  spur  gears  or  helical  gears,  an 
analysis  on  the  basis  of  plate  theory,  as  performed  by  Wellauer  and 
Seireg,33  is  appropriate.  However,  the  gears  used  in  the  program  dis¬ 
cussed  in  this  report  have  narrow  enough  face  widths  to  be  considered  as 
a  beam. 


LOAD  DISTRIBUTION 

An  area  of  design  interest  somewhat  similar  to  the  determination  of  dynamic 
factors  concerns  the  calculation  of  load  distribution  factors.  Because  of 
surface  asperities,  misalignment,  etc.,  the  tooth  load  may  not  be  evenly 
spread  across  the  tooth  face.  In  gear  design,  load  factors  are  used  to 
account  for  load  concentrations  on  the  tooth  which  otherwise  could  possibly 
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TOOTH  PAIR  COMPLIANCE,  in. /lb- -in. face  width(xl07) 


DEGREES  OF  GEAR  ROTATION  FROM  POSITION  WHERE 
A  TOOTH  PAIR  MESHES  AT  THE  PITCH  POINT 

Figure  8.  Tooth  Pair  Compliance  Versus  Contact  Position  as  Calculated 
by  Baud  and  Peterson. 
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MESH  COMPLIANCE,  in./lb.-in.  face  width(x  I07 ) 


20°-Full  Depth 
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DEGREES  OF  GEAR  ROTATION  FROM  POSITION  WHERE 
A  TOOTH  PAIR  MESHES  AT  THE  PITCH  POINT 


Figure  9.  Mesh  Compliance  Versus  Contact  Position  as  Calculated 
by  Baud  and  Peterson. 
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LOAD ,  In.  lb  f/ in.  face  width 


Figure  10.  Tooth  Deflection  Versus  Load  as 

Experimentally  Determined  by  Walker. 


DISTANCE  ON  LINE  OF  ACTION  FROM 
INTERFERENCE  POINT 

Figure  11.  Tooth  Deflection  Versus  Distance  on 

Line  of  Action  as  Determined  by  Walker. 
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cause  tooth  breakage.  Use  of  load  distribution  factors  is  appropriate  when 
considering  the  dynamic  response  of  geared  systems.  However,  in  the  narrow 
face-width  gears  used  in  this  study,  load  distribution  factors  were  found 
to  be  negligible. 

A  number  of  investigators,  e.g.,  Poritsky  et  al. j34  Weber,35 
Trbojevic,36*37  and  Conry,38  have  studied  load  distribution  factors. 
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GEAR  NOISE 

An  area  of  strong  interest  in  gear  dynamic  analysis  is  the  effort  to  pre¬ 
dict  and  control  gear  noise.  Most  of  the  work  performed  in  this  area  has 
been  on  a  semiempirical  basis.  Beuler  of  Volkswagen  presented  a  study  of 
geometrical  manufacturing  variables  on  gear  noise.39  Some  work  relating 
gear  noise  output  to  gear  quality  control  was  published  by  Bradley.40 

An  analytical  development  in  the  area  of  noise  was  presented  by  Welbourn.41 
In  it,  he  performed  a  Fourier  analysis  of  gear  tooth  errors  and  used  this 
information  to  determine  the  relative  effect  of  tooth  errors  on  gear  noise. 
Work  of  a  similar  nature  has  been  performed  by  Mechanical  Technology 
Incorporated.42"45  They  determined  Fourier  components  due  to  tooth  profile 
and  tooth  deflections  and  used  this  information  in  a  Holzer-type  torsion 
vibration  analysis  to  predict  vibration  and  noise. 


OTHER  WORK 

There  have  been  a  number  of  other  studies  which  cannot  be  categorized  in 
the  areas  discussed  previously.  Several  investigators  have  studied  tor¬ 
sional  dynamics  of  geared  systems.  Tordion  used  an  impedance  approach  to 
determine  natural  frequencies,  but  did  not  determine  instantaneous  loads.46 
Wang  developed  a  computation  procedure  to  predict  torsional  dynamic  re¬ 
sponse  which  accounted  for  all  system  components,  including  gears,  shafting, 
and  keyway  deflections.47  He  used  an  average  mesh  stiffness  value  in  his 
computations.  Fukuma,  Furukawa,  and  Aida  presented  similar  work,  which 
included  shaft  bending  deflections.48 

Another  interesting  study  was  presented  by  Wallace  and  Seireg,49  who  used  a 
finite  element  approach  to  study  the  dynamic  stress  and  deflection  in  a 
gear  tooth  due  to  a  moving  load  along  the  gear  tooth  profile. 


SUMMARY 


There  has  been  difficulty  in  the  development  of  gear  diagnostic  procedures 
because  of  a  lack  of  adequate  analytical  models.  This  deficiency  has  re¬ 
quired  gear  diagnostic  techniques  to  be  developed  almost  entirely  on  the 
basis  of  extensive  testing.  An  adequate  modelling  of  geared  system  dynamic 
response,  including  the  effect  of  tooth  surface  faults,  is  a  necessary  tool 
for  use  in  diagnostic  system  development.  While  portions  of  the  studies 
and  work  performed  in  the  area  of  gear  dynamics  may  be  applied  to  models  of 
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gears  with  tooth  surface  faults,  the  problem  as  a  whole  has  not  been 
adequately  analyzed*  Work  already  presented  may  be  used  to  gain  in¬ 
sight  and  direction  in  the  development  of  models  for  analysis  of  gears 
with  tooth  faults. 
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TESTING  APPARATUS 


Experimental  data  were  necessary  in  order  to  provide  a  method  of  verifying 
theoretical  predictions.  For  this  purpose,  a  test  stand  was  designed  and 
experiments  performed.  The  gearbox  used  on  the  test  stand  is  a  Generator 
Offset  Quill  Drive  from  a  UH-1  helicopter.  This  gearbox  is  instrumented 
to  provide  pertinent  data  from  tests  which  included  the  operation  of  both 
"good"  gears  and  gears  containing  faults  manufactured  on  the  tooth  surfaces. 


TEST  STAND 


Figure  12  shows  a  schematic  of  the  test  stand.  The  system  is  driven  by  a 
General  Electric  Model  26-G-440  Direct  Current  Electric  Dynamometer  which 
is  operated  as  a  motor.  This  motor  provides  up  to  20  hp  at  5000  rpm  and 
has  variable  speed  controls.  The  motor  is  coupled  to  a  Lebow  Model  1104 
Torquemeter  which  is  rated  at  500  in. -lb  and  9000  rpm.  The  torquemeter 
includes  a  built-in  magnetic  pickup  to  provide  a  speed  signal.  The  output 
shaft  of  the  torquemeter  is  then  coupled  to  the  input  spiral  bevel  gear  on 
the  gearbox.  Figure  13  shows  the  system  from  the  motor  to  the  gearbox 
(gearbox  on  the  right). 

The  gearbox  is  loaded  via  a  generator  and  resistor  bank.  The  generator 
used  is  a  Bendix  30B37-37-A  DC  Starter  Generator  which  is  rated  at  30  V, 

300  A  with  an  operation  range  of  4500-8000  rpm  and  is  the  same  one  used 
with  the  gearbox  in  general  helicopter  operation.  The  load  resistance  for 
the  generator  is  provided  by  twelve  1-ohm  resistors  wired  in  such  a  manner 
that  by  using  the  five  knife  switches  shown  on  the  control  board,  any  par¬ 
allel  combination  of  from  one  to  twelve  resistors  is  realizable.  Thus, 
load  resistance  is  variable  from  1  to  1/12  ohm.  Figure  l4  shows  the  con¬ 
figuration  of  the  generator  with  respect  to  the  gearbox,  and  Figure  15 
shows  the  load  resistance  bank,  A  rheostat  is  also  located  on  the  resistor 
bank  for  controlling  the  generator  field  and,  thus,  the  output  voltage. 


GEARS 


The  gearbox  contains  three  gears  as  shown  in  Figure  16.  The  driving  gear 
is  on  the  left, while  the  generator  driver  gear  is  on  the  right.  The  three 
gears  are  straight  spur  gears  having  the  design  parameters  shown  in  Table  I. 
The  overall  gear  ratio  through  the  gearbox  is  one;  i.e.,  the  input  speed  to 
the  gearbox  is  the  same  speed  as  that  of  the  generator  drive  shaft.  The 
ratio  between  each  gear  pair  is  1.35  and  0.746,  respectively.  The  design 
contact  ratio  is  1.75.  Gears  were  tested  both  in  "good"  condition  and 
with  faults  manufactured  into  them.  The  faults  were  designed  to  simulate 
pitch-line  pitting  and  were  produced  by  electrochemically  etching  a  line 
at  the  pitch  line  of  the  gear.  Figure  17  shows  a  sketch  of  the  faults. 

Two  sizes  of  faults  were  used:  a  0.020-inch-wide  fault  as  shown  in 
Figure  18  and  a  0.050-inch-wide  fault  as  shown  in  Figure  19.  A  gear 
which  had  a  0.020-inch-wide  fault  on  two  teeth  which  were  separated  by 
two  teeth  without  any  faults  was  also  tested.  Profile  checks  for  gears 
with  and  without  faults  are  shown  in  Appendix  I. 
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Figure  16.  View  of  Opene J 
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ith  0.020-Inch-Wide  Tooth  Fault 


TABLE  I.  GEAR  PARAMETERS 


Pinion 

Idler 

Gear 

Number  of  teeth 

4l 

55 

4l 

Diametral  pitch  (l/in.) 

10 

10 

10 

Pressure  angle  (deg) 

20 

20 

20 

Pitch  diameter  (in.) 

4.10 

5-5 

4.10 

Face  width  (in.) 

0.20 

0.40 

0.20 

Circular  pitch  (in.) 

0.3142 

0.3142 

0.3142 

TRANSDUCERS 


A  number  of  parameters  were  monitored  during  the  test  program.  A  wide 
variety  of  accelerometer  positions  are  possible  on  the  gearbox.  However, 
three  positions  appeared  to  provide  the  greatest  amount  of  information  ,  and 
these  were  used  more  extensively  than  others.  Figure  20  shows  a  Brilel  and 
Kjaer  4344  accelerometer  which  was  stud-mounted  near  the  input  bearing  of 
the  gearbox.  Figure  21  shown  an  Endevco  2215  accelerometer  which  was 
epoxied  to  the  face  of  the  gearbox  housing.  Figure  22  shows  a  Bruel  and 
Kjaer  4344  accelerometer  which  was  stud-mounted  near  the  output  bearings  of 
the  gearbox. 

Strain  gages  were  mounted  at  the  root  fillets  of  specific  gear  teeth  to  pro¬ 
vide  information  on  the  dynamic  strain  occurring  during  meshing.  Figure  23 
shows  a  typical  mounting  configuration  of  three  strain  gage  bridges 
on  a  gear.  Active  gages  were  mounted  on  the  dedendum  of  three  consecutive 
gears.  Dummy  gages  were  mounted  in  the  middle  of  teeth  away  from  the  active 
gages  as  well  as  on  the  opposite  face  of  the  gear  (not  shown  in  Figure  23). 
The  strain  gage  signals  are  transmitted  from  the  gear  via  a  set  of  18 
channel  Superior  Carbon  Sliprings.  The  position  of  the  slipring  assembly  is 
shown  in  Figure  14. 

Torque  and  speed  were  monitored  with  signals  provided  by  the  torquemeter. 
Output  power  was  monitored  with  a  voltmeter  and  the  knowledge  of  load  re¬ 
sistance. 


SIGNAL  CONDITIONING  AND  ANALYSIS  EQUIPMENT 


Figure  24  shows  some  of  the  signal  conditioning  and  analysis  equipment  used 
in  testing.  Kistler,  Kristal,  and  Unholtz  Dickie  charge  amplifiers  are  used 
with  the  accelerometer  inputs.  Honeywell  Accudata  strain  gage  signal  con¬ 
ditioners  and  an  input  amplifier  are  used  with  t.ie  strain  and  torque  signals. 
An  EPUT  counter  is  used  with  the  tachometer  signal. 
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Mounting  Position  of  Accelerometer  Near  Input  Bearing 


Data  were  either  analyzed  on-line  or  tape-recorded  and  analyzed  later.  A 
Precision  Instruments  6100  FM  tape  recorder  with  three  available  channels 
was  used  for  the  tape  recording. 

Much  of  the  data  analysis  was  done  with  a  Federal  Scientific  UA-500 
Ubiquitous  Spectrum  Analyzer.  This  instrument  provides  the  Fourier  com¬ 
ponents  of  signals  fed  into  it  either  directly  from  the  instrumentation 
during  test  runs  or  from  taped  data. 

Transfer  function  and  impedance  tests  were  performed  on  the  test  stand  with 
the  use  of  the  Spectral  Dynamics  Model  SD1002E  Mechanical  Impedance  Systems, 
shown  in  Figure  25.  Excitation  was  provided  with  an  MB  Electronics  50-pound 
shaker. 
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THEORETICAL  ANALYSIS 


A  theoretical  analysis  of  the  geared  system  which  was  being  used  for  the 
generation  of  test  data  was  performed.  This  analysis  resulted  in  a  model 
which  provides  both  frequency-domain  and  time-domain  response  information 
for  the  system  when  all  the  gears  are  in  good  condition  and  also  when  gears 
contain  tooth  surface  faults.  The  model  consist j  of  a  torsional  dynamic 
analysis  of  the  geared  system  and  a  separate  determination  of  the  gear  mesh 
stiffness  functions.  In  both  analyses,  modifications  are  incorporated  to 
account  for  the  effects  of  tooth  surface  faults. 


TORSIONAL  DYNAMIC  ANALYSIS 


The  gears  contained  in  the  gearbox  are  modeled  as  shown  in  Figure  26,  with 
the  equivalent  rotational  inertias  of  the  gears  coupled  by  a  linear  spring 
connecting  their  base  circles.  The  spring  stiffness  is  a  function  of  the 
meshing  position  of  the  gear  teeth. 


Figure  26.  Torsional  Model  of  a 
Gear  Pair. 
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To  determine  the  overall  response  of  the  geared  system,  a  lumped  parameter 
torsional  analysis  of  the  system  shown  in  Figure  27  is  performed.  The 
following  are  the  pertinent  parameters  shown  in  the  figure  and  used  in  this 
analysis: 

Ji  inertia  of  shafting  and  input  gear,  in.-lb»s2 

J2  inertia  of  spur  gear  pinion,  in.-lb*s2 

J3  inertia  of  idler  spur  gear,  in.-lb«s2 

J4  inertia  of  spur  gear,  in.-lb*s2 

•T5  inertia  of  generator,  in.-lb»s2 

K33  linear  spring  coefficient  of  pinion- idler  mesh, 

lb/in. 

B23  linear  damping  coefficient  of  pinion-idler  mesh, 

lb/in.  *s_1 

K34  linear  spring  coefficient  of  idler-gear  mesh,  lb/in. 

B34  linear  damping  coefficient  of  idler-gear  mesh, 

lb/in.*s_1 

Ktia  torsional  spring  coefficient  of  input  shafting, 

in. /lb*  rad-1 

B.l2  torsional  damping  coefficient  of  input  shafting, 

in.-lb/rad*s-1 


K+45  torsional  spring  coefficient  of  spur  gear  shaft, 

in. -lb /rad 

B.  45  torsional  damping  coefficient  of  spur  gear  shaft, 

in.-lb/rad's"1 

9i  >  02,  03,  84,  05  angular  rotation  of  inertias,  rad 
eg  generator  output  voltage,  V 

ig  generator  outpit  current,  A 

Rg  generator  output  load,  fi 

K_  generator  voltage  constant,  v/rad^s"1 

D 

generator  current  constant,  in.-lb/A 
r-  radius  of  pinion  base  circle,  in. 
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'MS/i.'  4 


Ktl2’ 


torsional  Model  of  the  Geared  System. 
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V. 


The  analysis  of  the  torsional  mode  of  vibration  of  the  system  was  performed 
since  this  was  expected  to  be  the  primary  mode  of  vibration  encountered. 

To  translate  from  the  torsional  response  of  the  system  to  lateral  accelera¬ 
tions  monitored  on  the  gearbox,  experimentally  generated  transfer  function 
information  between  the  gear  mesh  and  the  housing  is  used.  The  system  is 
excited  by  both  the  input  driving  torque  and  the  variations  of  the  mesh 
stiffness  which  occur  as  a  function  of  gear  rotational  position. 

The  electric  generator  is  considered  to  operate  with  a  constant  field  volt¬ 
age  and  negligible  armature  inductance.  Load  resistance  is  kept  at  a  con¬ 
stant  value. 

Torque  balances  for  each  of  the  lumps  are  written  as  follows: 

Jl  ©i  *  Ktl2(®3“®l)  +  ®tl2(^2“®i)  +  (5) 

•  •  •  • 

J9O9  *  -K-^i  9  (69-61  )  -  Btl2(^»“^l  )  +  (R3@3-R209  )  +  B23R?  ) 

(6) 

J303  0  “K03R3  (R30'3-R3$?  )  “  B93R3 ( R3 03“R9 0? )  +  K34R3 ( R404— R3 63  ) 

+B34ft3(B4^4"B3^3)  (7) 

J4O4  =  — K34R4 ( R4O4— R3 $3 )  +  B34R3  (R404-R30s)  +  ^45(05-^4)  +  6^43(03-04) 

(8) 


J 5^5  a  ”K^45(^3"^4)  "  8^.45(05-04)  -  Ktpig  (9) 

Use  of  Kirchoff's  Law  on  the  generator  loop  yields 

Kg 05  -  igPg  =0  (10) 

It  is  convenient  to  first  rearrange  the  equations  in  the  following  form: 

01  =  r  K^  2  ( 03 -01  )  +  B.^  2(0?“0i  )  +  TX1  (ll) 

Ji  L 

0?  =  y—  ^K?3Rf  ^  ^.03-02^  +  B93R3  ^  ^*^3”^?^  "  K^12(02“0i)  “  B^12(®2”^l)j 

(12) 

03  =  j-  ^K34R3R4  ^  ©4"  ^03  ^  +  B34R3R4  ^  04-  ^303  ^  "  K23R3R9  ^  j^l03“02^ 

-B23R3R?  ^  ^03"0?^  ^  (13) 
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O5-O4 )  “  K34R| 

(*•-  kk)} 


(04-  §^3) 


-B34R4 


(14) 


j-  £-^45(05-04)  -  B^45(05-04)  -  KijiigJ 


ig  -  i-Kge5 


(15) 

(16) 


Parameters  representing  a  relative  rotation  between  the  system  inertias  are 
now  defined  as 


*  03  -  Qi 

(17) 

:  “  0? 

(18) 

1  04  —  — ^03 

(19) 

r4 

■05-04 

(20) 

The  utility  of  these  parameters  is  seen  if  we  consider  the  response  of  the 
system  rotational  angles.  Figure  28  shows  a  typical  response  curve  of  one 
of  the  angles ,  0X .  The  response  is  a  perturbation  about  a  steadily  in¬ 
creasing  ramp  function.  This  ramp  increases  at  a  rate  corresponding  to  the 
input  speed.  The  perturbations  are  the  dynamic  effects  of  the  system's 
stiffnesses,  inertias,  and  damping.  Unless  the  driving  shafts  reverse 
direction,  the  rotational  angles  will  continue  to  grow.  As  time  progresses, 
the  percentage  variation  due  to  the  dynamic  perturbations  will  become 
smaller  and  smaller.  Computationally,  the  perturbations  become  more  diffi¬ 
cult  to  resolve.  And,  since  vibrations  caused  by  this  relative  motion 
between  the  inertias  are  the  desired  results,  it  is  these  parameters  which 
are  of  primary  interest. 

The  new  parameters  may  be  differentiated  once  with  the  result 


1^1  =  0?  -  01 

(21) 

\Jfj_  »  ~a,03  -  0? 

R-» 

(22) 
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*3  »  04  -  Sao3 

K4 


(23) 


1r4  =  Os  ~  64 

Combining  Equations  (11-24)  yields 

=  j-  (  2\|rT  +  +  T,  ) 


(24) 


(25) 


^2  =  (^23^2^2  +  B23R2t2  “  ^tl  2'lT'L  ”  Bti2^i ) 

5a@3  =  ^  —  (  K34R3R4t3  +  B34R3R4\if3  -  K23R3R2\|r2  -  B?3R3R3v|r3  ) 
h2  r2  J3 

04  =  (  ^45^4  +  B-t45\if4  -  K34R4\(f3  -  B34R4\|f3) 

03  =  (  -  K-t45\Jf4  -  B‘t43'!/4  -  K>pig  ) 


ig  - 


In  the  system  analyzed, 


Rj»  =  R4 


Thus, 

R3  =  Rg 

R5  R4 

Using  Equation  (31)  and  differentiating  Equations  (20-23)  we  find: 

=  ^5  ” 


(26) 

(27) 

(28) 

(16) 

(30) 

(31) 

(32) 


b  =  — °3  -  b  (33) 

r2 

b  =  o4  -  |a.03  (34) 

t4  =  05  "  04  (35) 

Equations  (24-29)  may  now  be  combined  with  Equations  (32-35).  It  should  be 
noted  that  these  equations  contain  position-varying  coefficients  due  to  the 
gear  stiffness  function.  This  fact  eliminates  the  possibility  of  an  easily 
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determined  closed-form  solution  for  system  response.  Therefore,  the  solu¬ 
tion  is  obtained  by  using  a  numerical  integration  scheme  via  CSMP  (Con¬ 
tinuous  System  Modelling  Program)  on  the  IBM  370/165,  as  explained  in 
Appendix  II. 


GEAR  STIFFNESS  CALCULATIONS 

As  previously  discussed,  the  mesh  stiffness  of  a  gear  varies  as  a  function 
of  gear  rotational  position.  This  variation  is  a  result  of  both  the  complex 
gear  tooth  shape  and  the  load-sharing  which  is  exhibited  between  gear  tooth 
pairs.  The  approach  which  was  taken  to  determine  the  too+h  stiffness  func¬ 
tion  was  a  combination  of  an  elastic  strain  energy  approach  and  a  classical 
Hertzian  contact  procedure  developed  by  Weber32  and  used  by  Richardson, 19 
Wang,47  and  others.  The  analysis  considers  all  tooth  stiffness  to  be  con¬ 
centrated  in  the  tooth  mesh;  i.e.,  the  gear  web  stiffness  is  considered 
negligible. 

The  generalized  approach  is  to  separately  determine  the  tooth-bending  and 
Hertzian  compliances.  These  are  then  summed  to  provide  total  tooth  pair 
compliance.  Load-sharing  considerations  are  then  used  to  determine  the 
total  mesh  stiffness  function. 

In  calculating  the  stiffness  functions,  a  static  analysis  is  performed  with 
respect  to  tooth  geometry.  The  teeth  are  assumed  to  mesh  solely  along  the 
line  of  action.  Any  dynamic  actions  or  tooth  errors  which  will  cause  motion 
off  the  line  of  action  are  not  included  in  this  analysis.  Such  meshing 
action  would  cause  changes  in  load-sharing  and  would  require  an  iterative 
solution  procedure  in  order  to  locate  the  actual  points  of  tooth  contact. 


Bending  Compliance 

Figure  29  shows  the  tooth  model  as  analyzed.  The  force  FL  is  imposed  on 
the  tooth  at  a  given  contact  position.  Based  on  the  contact  position,  an 
angle  0  (the  angle  between  the  direction  of  force  application  and  a  perpen¬ 
dicular  to  the  tooth  centerline)  and  the  height  YL  above  the  tooth  root 
circle  are  determined.  Using  these  two  parameters,  the  load  force  is  re¬ 
solved  into  a  moment  M,  a  shear  force  Q,  and  a  normal  force  N,  all  refer¬ 
enced  to  the  tooth  centerline.  These  are  expressed  as  follows: 


Coordinates  used  are  Y 
position  on  the  tooth) 
point) . 


Q  =  Fl  cosfl  (36) 

N  -  Fl  sin$  (37) 

M  =  Q(Yl  -  Y)  (38) 

(the  distance  from  the  tooth  root  circle  at  any  given 
and  X(Y)  (the  half -thickness  of  the  tooth  at  any 
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The  deflection  work  is  equated  to  the  strain  energy  to  give 

fF^B  -!/fl“tl/,|g!4Y+  l/S41 


where 


(39) 


0 

6B  =  tooth  deflection  due  to  bending,  in. 

Fl  «  force  imposed  upon  the  tooth,  lb 
M  »  moment  on  tooth  resolved  to  tooth  centerline,  in. -lb 
E  »  Young's  modulus,  lb/in,2 
I  a  tooth  moment  of  inertia,  in.4 

k  a  correction  factor  used  to  account  for  nonuniform  shear 
distribution 

Q  =  shear  on  tooth  resolved  to  tooth  centerline,  lb 
A  a  tooth  cross-sectional  area,  in.2 

N  =  normal  force  on  tooth  referenced  to  tooth  centerline,  lb 


I  - 


W(2X): 


12 

A  =  W(2X) 

where  2X  =  thickness  of  tooth,  in. 

Combining  Equations  (38-4l)  and  dividing  by  1/2  FL  gives 

6b  =  k /  dY  +  k  I  dY  +  k  I  dY 


(40) 

(41) 


(42) 


A  new  parameter  FQ  is  defined  to  represent  the  normalized  tooth  force 

F 

0  "  W 


(43) 


where  W  =  facewidth,  in. 
Combining  Equations  (36-42)  gives 


yl  yl  yl 

_  12F0  cos2(g)  f  (YL-Y)2  kF0  cos2 (9)  f  _dy_  F0  sin2 (9)  r  _dy 
3  '  ^  E  l  (2X)3  G  I  (2X)  E  J0  (2X 


(2X) 

(44) 
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A  shear  correction  factor  of  1.2  (that  used  for  rectangular  beams)50  is 
also  used  with  the  resulting  bending  compliance  expression 


12  cos2(0)  (Yt-Y)2  T  1.2  cos2(0)  sin2(9)  dY 

-  I  — i - dY  +  - + -  I  - 

E  l  (2X) 3  G  E  J0  (2X) 

L  (45) 


It  should  be  noted  that  this  compliance  is  defined  with  respect  to  a  given 
tooth  meshing  position  through  the  parameters  Y,  X,  YL>  and  9. 


Hertzian  Compliance 

The  equations  used  to  describe  the  Hertzian  or  contact  deflection  and, 
thus,  Hertzian  compliance  are  a  modification  of  those  used  to  describe 
contact  between  two  cylinders.51*52  First,  the  width  of  the  contact  strip 
(b)  between  the  two  contacting  teeth  is  found: 


where  b  =  width  of  contact  strip,  in. 

rG  =  distance  from  contact  point  to  tangent  line  of  base  circles 
on  gear,  in. 

rp  =  distance  from  contact  point  to  tangent  line  of  base  circles 
on  pinion,  in. 
v  =  Poisson's  ratio 

The  Hertzian  deflection  is  now  expressed  as 


6 


2F0(l-v2)  f  2hG  2hp  v  1 

-  fin - +  fin -  -  -  J 

rtE  L  b  b  1-vJ 


(47) 


where  hG  =  distance  from  contact  point  to  the  tooth  centerline  of 
gear  along  the  angle  of  tooth  loading,  in, 
hp  =  distance  from  contact  point  to  the  tooth  centerline  of 
P  pinion  along  the  angle  of  tooth  loading,  in. 


Expressing  this  in  terms  of  compliance, 


C 


H 


2  ( 1-v ) 7  [  2hr 

— —  t  + 


Or 


2h_,  v 

“F"  "  T7v. 


(48) 


56 


It  should  be  noted  in  the  expressions  for  Hertzian  compliance  that  a  non¬ 
linearity  with  respect  to  load  force  is  involved;  i.e.,  the  width  of  the 
contact  strip  is  a  function  of  load  force.  In  calculation  of  this  campli 
ance,  the  mean  value  of  tooth  force  is  used  to  find  the  contact  width. 

Mesh  Stiffness 


Equations  (45-47)  are  solved  to  provide  compliance  values  as  a  function  of 
meshing  position.  Details  of  geometric  relationships  and  computations  are 
given  in  Appendixes  III  and  IV,  To  find  overall  mesh  stiffness,  an 
expression  is  first  used  for  gear  pair  total  compliance;  i.e., 

CTotal  =  CB  gear  +  CB  pinion  +  CH  (49) 

The  gear  pair  stiffness  is  the  inverse  of  this  compliance;  i.e., 


^pair 


1 

cTotal 


(50) 


To  determine  the  complete  mesh  stiffness  function,  tooth  pair  load  sharing 
must  be  introduced.  Load  sharing  occurs  at  the  distance  of  the  normal  pitch 
from  either  end  of  the  line  of  action  as  shown  in  Figure  30.  This  allows 
determination  of  the  positions  at  which  more  than  one  tooth  pair  are  in  con¬ 
tact  and,  thus,  where  load-shoring  must  be  included  in  the  analysis. 


Figure  30»  Tooth  Load  Sharing  Positions, 
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When  more  than  one  tooth  pair  are  in  contact,  the  constraint  that  the  rigid 
body  motion  of  each  pair  must  be  equal  is  used. 


&  =  =  &2 


FTctal  =  Fi  +  F? 


Ft 

h 
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FTotal 


force  on  first  tooth  pair,  lb 
stiffness  of  first  tooth  pair,  lb/in. 
deflection  of  first  tooth  pair,  in. 
force  on  second  tooth  pair,  lb 
stiffness  of  second  tooth  pair,  lb/in. 
deflection  of  second  tooth  pair,  in. 
total  force  on  mesh,  lb 


Combining  Equations  (51-54)  ,  we  get 


FTotal 

6 


Kx  +  Kg 


(51) 

(52) 

(53) 

(54) 


(55) 


Thus,  to  compute  the  effective  mesh  stiffness,  the  individual  stiffnesses 
are  added  in  the  manner  of  springs  in  parallel.  Figure  31  is  a  graphical 
representation  of  this  procedure.  This  stiffness,  which  is  used  in  the 
torsional  dynamic  analysis,  it  represented  as  a  periodic  function  of  gear 
rotation. 


EFFECT  OF  FAULTS 


The  previous  analysis  referred  to  gear  in  "good"  condition,  i.e.,  exhibiting 
no  tooth  surface  faults.  To  predict  the  response  of  gears  having  faults, 
two  modifications  were  made  in  the  analysis.  A  recalculation  of  stiffness 
to  include  the  fault  effect  on  the  stiffness  functions  was  m  ie.  As  an 
approximation  to  the  tooth  impact  which  occurs  at  the  mesh,  a  torque  impact 
was  also  introduced  into  the  torsional  dynamic  equations. 


Stiffness  for  Tooth  With  Fault 


The  first  manifestation  of  the  tooth  surface  fault  is  introduced  in  terms 
of  the  tooth  thickness  (2X).  The  thickness  function  is  modified  by  the 
amount  of  material  which  is  missing  due  to  the  fault  (see  Figure  32).  This 
change  in  thickness  causes  a  change  in  the  gear  tooth  bending  compliance. 
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Figure  32.  Model  of  Tooth  Fault. 
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The  Hertzian  compliance  analysis  is  also  affected  by  the  modified  tooth  pro- 
file.  Whenever  the  fault  is  in  mesh,  the  width  of  the  contact  strip  is 
modified  by  the  amount  of  overlap  between  the  fault  and  the  original  con¬ 
tact  strip.  If  the  overlap  of  the  fault  is  greater  than  the  width  of  the 
contact  strip,  it  is  also  assumed  that  contact  is  momentarily  lost.  When 
this  effect  occurs,  the  stiffness  function  is  set  equal  to  zero.  The  three 
considerations  discussed  as  effects  of  a  fault  are  incorporated  into  the 
calculation  of  a  second  stiffness  function  for  a  tooth  with  a  fault  (see 
Appendix  IV).  Whenever  contact  is  lost  in  the  mesh,  it  is  assumed  that 
damping  effects  are  also  lost.  Thus,  when  mesh  stiffness  goes  to  zero, 
mesh  damping  also  goes  to  zero  in  the  torsional  analysis. 


Torque  Impact 

Whenever  contact  is  lost  because  of  a  fault,  the  gears  can  freewheel;  i.e., 
the  pinion  tends  to  accelerate  since  it  has  no  resisting  torque,  while  the 
idler  gear  tends  to  decelerate  because  of  the  loss  of  its  driving  torque. 
When  contact  is  reestablished,  an  instantaneous  velocity  difference  between 
the  two  teeth  is  present  which  causes  an  impact  at  the  mesh. 

This  impact  is  approximated  as  a  torque  upon  the  gear  mesh.  The  change  in 
angular  momentum  between  the  gears  determines  the  impact  amplitude. 


rIMP 


JaA^»(t*)  -  JVA,(t*) 
tIMP 


(56) 


where 


TIMP 

A0?(t*) 

A^(t*) 

TIMP 


amplitude  of  torque  Impact,  in. -lb 

change  in  velocity  of  pinion  at  time  of  contact,  rad/s 

change  in  velocity  of  gear  at  time  of  contact;  rad/s 
length  of  time  of  impact,  s 


The  impact  is  entered  in  the  analysis  as  torques  caused  by  equal  and  oppo¬ 
site  forces  acting  on  the  sprirg  connecting  the  two  inertias  of  the  gear 
mesh  ( Jg > J3 ) ,  These  torques  are  included  in  Equations  (6)  and  (7)  whenever 
tooth  contact  reaches  the  end  of  the  fault  width. 

The  impact  is  expected  to  assume  a  triangular  shape  with  respect  to  time 
(see  Figure  33).  However,  because  of  the  fixed-interval  integration  rou¬ 
tine  used,  the  rectangular  model  of  the  impact,  as  shown  by  the  dashed 
lines,  is  used.  The  energy  content  of  the  impact  is  unchanged  with  this 
model  and  more  repeatability  of  results  is  obtainable.  This  is  the  result 
of  the  fixsd-interval  integration  procedure  used  in  the  simulation  which 
does  not  necessarily  perform  the  computation  at  the  same  point  on  the 
impact  each  time  the  fault  enters  the  mesh. 
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SUMMARY 
■■■  ■ 

s? 

A  set  of  torsional  dynamic  equations  was  written  to  describe  relative 
motion  between  the  inertias  of  the  model.  The  stiffness  function  of  the 
gear  mesh  was  also  described  via  a  strain  energy  analysis.  Modifications 
to  the  stiffness  function  and  the  inclusion  of  a  torque  impact  were 
described  to  account  for  gear  tooth  faults  in  mesh, 
i 
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EXPERIMENTAL  RESULTS 


Experiments  were  performed  using  the  equipment  described  in  the  section  on 
test  apparatus.  Acceleration,  torque,  speed,  and  strain  data  were  monitored 
for  a  variety  of  test  conditions  and  a  number  of  different  gears.  The 
standard  test  procedure  used  an  input  speed  of  4500  ±  10  rpm,  with  an  out¬ 
put  power  level  of  8  ±  0.5  kW.  These  conditions  fall  within  the  operating 
range  of  the  gearbox  while  in  normal  helicopter  usage.  Tests  were  also 
performed  at  other  speed  and  load  conditions  to  determine  these  effects  on 
overall  response  data. 


GEARS 


A  total  of  seven  gears  (spur  pinion,  generator  drive,  transmission)  were 
used  in  the  test  program.  These  gears  are  manufactured  to  be  symmetric 
and  reversible,  thus  giving  the  possibility  of  14  distinct  sets  of  tooth 
surfaces.  Table  II  lists  the  gears  by  serial  number  and  testing  condition. 
The  notations  R  and  L  refer  to  the  right  and  left  sides  of  the  gear  tooth, 
respectively  (viewed  with  the  serial  number  facing  up) .  The  table  describes 
whether  a  given  gear  side  was  run,  what  the  surface  conditions  of  the  teeth 
were  (i.e.,  if  the  gear  contained  a  fault  or  not),  and  if  the  gear  had 
strain  gages  mounted  on  it. 


TABLE  II.  GEARS  TESTED 


Gear  No. 

Condition 

B12  - 

1659 

R 

Run,  Good 

L 

Run,  Small  Fault 

W12  - 

3231 

R 

Not  Run 

L 

Run,  Large  Fault, 

Gaged 

B12  - 

1475 

R 

Not  Run 

L 

Run,  Good,  Gaged 

W12  - 

3204 

R 

Run,  Good 

L 

Run,  Double  Fault 

W12  - 

3191 

R 

Run,  Good,  Gaged 

L 

Not  Run 

W12  - 

3183 

R 

Not  Run 

L 

Run,  Good 

W12  - 

3115 

R 

Run,  Small  Fault, 

Gaged 

L 

Run,  Good 

TIME- DOMAIN  ACCELERATION  AND  TORQUE  RESPONSE 


Gearbox  acceleration  signals  are  among  the  parameters  of  most  interest, 
both  with  respect  to  experimental  information  and  comparison  to  theoretical 
results.  Acceleration  data  presented  was  monitored  at  one  of  three  posi¬ 
tions  as  shown  in  Figures  20,  21,  and  22  (near  the  gearbox  input  bearing, 
on  che  face  of  the  gearbox,  and  near  the  output  bearing).  The  position 
winch  provided  the  greatest  response  to  faults  was  at  the  input  bearing 
3ince  this  was  the  closest  monitoring  position  to  the  gear  mesh  which  had 
faulty  gears  implanted  in  it. 

Figure  3^  shows  typical  torque  and  acceleration  signals  for  the  system 
when  all  the  gears  were  in  good  condition.  The  top  trace  in  the  figure  is 
the  torque  signal.  The  straight  line  is  the  zero  torque  level.  The  input 
torque  has  a  mean  level  with  a  superimposed  oscillation  at  the  input  shaft 
frequency.  This  oscillation  was  deduced  to  be  the  result  of  mounting  pro¬ 
cedures  and  fixtures  used  in  the  test  stand.  Accordingly,  this  information 
was  used  as  an  input  in  the  system  simulation.  The  bottom  trace  in  the 
figure  is  for  the  input  bearing  accelerometer,  the  top  trace  in  figure  (b) 
is  the  gearbox  housing  accelerometer  signal,  and  the  bottom  trace  in  figure 
(b)  is  the  output  bearing  accelerometer  signal.  The  primary  oscillatory 
component  of  these  signals  occurs  at  the  gear  meshing  frequency.  Some 
"beating"  or  modulation  effects  at  the  shaft  rotational  frequency  may  also 
be  noted. 

Figure  35  contains  data  similar  to  that  of  Figure  34  for  a  gearbox  contain¬ 
ing  a  gear  with  a  small  ( 0.020-inch -wide)  tooth  surface  fault.  The  top 
trace  in  each  of  the  scope  photos  is  the  tachometer  signal.  Some  small 
spikes  occurring  at  a  repetition  rate  which  corresponds  to  the  shaft  rota¬ 
tional  frequency  may  be  seen  in  the  acceleration  signal.  The  repetition 
rate  corresponds  to  the  frequency  at  which  the  fault  enters  mesh.  The 
presence  of  the  spike  is  evident  but  not  very  large.  Figure  35  presents 
data  for  the  gearbox  when  a  gear  with  a  large  (~0, 05 -inch -wide)  fault  had 
been  implanted.  The  presence  of  this  spike  is  more  evident. 

Differences  between  good  gears  and  those  with  faults  are  seen  in  the  accel¬ 
eration  data.  However,  because  of  the  magnitudes  of  the  spikes,  it  is 
questionable  as  to  whether  the  presence  of  these  spikes  could  be  used  solely 
as  a  gear  condition  indicator  or  whether  they  would  be  dismissed  as  anomalies 
in  the  test  procedure  if  the  gear  condition  was  unknown.  The  use  of  time 
summation  as  an  analysis  procedure  has  been  suggested  for  use  in  such  cases.8 
However,  based  on  the  predominance  of  energy  in  the  signal  at  the  gear  mesh¬ 
ing  frequency  and  its  harmonics,  it  is  quite  doubtful  whether  a  summation 
at  the  shaft  period  would  produce  useful  information. 


ACCELERATION  FREQUENCY  SPECTRA 

The  major  portion  of  the  analysis  of  data  was  performed  using  frequency 
spectra  of  the  input  bearing  acceleration.  These  frequency  analyses  were 
obtained  with  the  use  of  a  real-time  analyzer  for  data  taken  both  directly 
from  the  test  stand  and  from  taped  data  which  were  analyzed  off-line. 
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Figure  36.  Time  Response  of  Acceleration  and  Torque 
(Large  Fault). 
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Figure  37  shows  typical  spectra  for  the  input  bearing  accelerometer  when 
the  gearbox  contained  all  good  gears.  The  input  shaft  frequency  is  75  Hz 
and  the  fundamental  gear  meshing  frequency  is  3075  Hz,  In  Figure  37a 
(0-10  KHz),  the  fundamental  component  and  the  first  two  harmonics  may  be 
seen.  Also  evident  are  sidebands  (primarily  about  the  first  and  second 
mesh  harmonics).  It  is  quite  evident  that  the  primary  composition  of  the 
signal  is  due  to  the  gear  meshing  frequency  and  harmonics.  Figure  37b 
(0-50  kHz)  showo  more  of  the  harmonic  train.  It  is  interesting  to  note 
that  the  first  harmonic  has  a  larger  amplitude  than  the  fundamental  gear 
meshing  frequency,  something  which  is  contrary  to  the  expected  result. 

Figure  38  shows  spectra  when  the  gearbox  contained  a  gear  having  a  small 
(0.020 -inch -wide)  fault.  With  a  fault  present  the  amplitude  at  the  mesh¬ 
ing  frequency  is  increased.  This  amplitude  is  now  also  greater  than  that 
of  the  first  harmonic.  There  is  the  appearance  of  a  larger  amount  of 
energy  at  multiples  of  the  shaft  frequency  (the  sawtooth  effect  which  may 
be  observed).  Figure  39  shows  spectra  when  a  gear  containing  a  large 
fault  (0.050  inch  wide)  was  tested.  The  same  effects  as  noted  for  the 
small  fault  may  be  observed,  with  the  sawtooth  effect  having  still  higher 
amplitudes.  Figure  40  presents  spectra  obtained  when  the  gear  tested  con¬ 
tained  two  faulty  teeth.  Again,  the  increases  in  amplitude  and  sawtooth 
effect  are  evident.  Gears  with  faults  on  more  teeth,  as  might  be  indica¬ 
tive  of  normal  gear  tooth  failures,  were  not  tested  because  of  the  magni¬ 
tude  of  machining  required. 

Another  effect  which  is  especially  evident  in  Figure  39  is  a  periodicity 
in  the  sawtooth  effect.  An  enveloping  periodic  amplitude  change  occurs 
with  a  spacing  which  corresponds  to  the  lowest  torsional  resonance  of  the 
system.  This  is  particularly  noticeable  at  frequencies  between  the  gear 
meshing  frequency  and  its  first  harmonic  in  Figure  39. 


Figures  4l-44  display  the  variation  in  test  to  test  amplitudes  at  various 
selected  frequencies.  Amplitudes  at  the  meshing  frequency  and  its  harmonics 
are  displayed  in  terms  of  decibels  referenced  to  0.05  g  (0.05  g  =  0  dB) . 

The  vertical  lines  on  the  plots  of  the  amplitude  of  meshing  frequency  and 
its  first  two  harmonics  are  mean  values  for  good  and  faulty  gears,  respec¬ 
tively.  The  mean  value  of  the  meshing  frequency  is  approximately  9  dB 
higher  for  faulty  gears  than  for  good  gears.  Table  III  displays  the  am¬ 
plitudes  of  the  mean  values  and  also  the  mean  amplitude  plus  and  minus  one 

TABLE  III.  VARIATIONS  IN  SPECTRAL  AMPLITUDES  (dB) 

OF  GEARS  TESTED 


Good  Gears_  _  Faulty  Gears 


Frequency 

X 

X  -  CT 

X  +  o- 

X 

X  -  CT 

X  +  o- 

Gear  Meshing 

35.3 

31.9 

37.7 

44.5 

42.2 

46.2 

First  Harmonic 

42.5 

32.1 

47.1 

39.4 

34.8 

42.4 

Second  Harmonic 

26.8 

20.2 

30.5 

27.6 

21.4 

31.2 
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FREQUENCY  (kHz) 


Figure  38.  Acceleration  Spectra  at  Input  Bearing 
(Small  Fault) . 


ACCELERATION  (g) 


FREQUENCY  (kHz) 


Figure  40.  Acceleration  Spectra  at  Input  Bearing 
(Double  Fault). 
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Figure  42.  Distribution  of  Spectral  Amplitudes  at  First  Harmonic  and  Sidebands  at  Input 
Bearing  Sensor  Location. 
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Figure  44.  Distribution  of  Spectral  Amplitudes  at  Third  Through  Seventh  Harmonic  at  Input 
Bearing  Sensor  Location. 


standard  deviation.  The  changes  caused  by  gear  faults  are  not  very  distinct 
in  the  harmonics  of  the  meshing  frequency  or  sidebands  as  data  points  tend 
to  overlap.  While  the  presence  of  a  tooth  surface  fault  undoubtedly  affects 
spectral  amplitudes  of  these  frequencies,  other  operating  parameters  of  the 
gearbox  must  also  strongly  affect  them.  This  suggests  that  the  primary 
information  with  respect  to  fault  indication  lies  in  the  amplitude  of  the 
gear-meshing  frequency  and  the  presence  of  energy  at  multiples  of  the  rota¬ 
tional  frequency,  the  sawtooth  effect. 

Figure  45  shows  typical  spectra  for  good  gears  when  load  was  set  to  70 %  of 
the  standard  operating  load.  No  great  amount  of  deviation  from  the  effects 
noted  under  full  load  were  observed.  Figures  46  ,  47,  and  48  show  the  varia¬ 
tions  at  various  spectral  frequencies  for  good  and  bad  gears  at  the  reduced 
load  condition.  Again  amplitudes  at  the  meshing  frequency  and  first  har¬ 
monics  appear  higher  for  faulty  gears,  but  no  clear-cut  pattern  is  evident 
in  the  sideband  information.  Other  data  were  taken  at  40$  of  the  standard 
operating  load.  The  spectra  contained  an  increase  of  harmonic  content,  and 
gears  were  much  noisier  at  this  operating  condition. 

Figures  49  and  50,  respectively,  display  spectra  of  acceleration  monitored 
at  the  generator  output  bearing  when  the  gearbox  contained  all  good  gearr 
and  when  a  faulty  gear  was  implanted.  Differences  may  be  noted  in  the 
amplitudes  at  the  meshing  frequency  and  harmonics  when  a  faulty  gear  was 
present.  The  faulty  gear  also  shows  some  evidence  of  the  sawtooth  effect. 
Figures  51,  52,  and  53  show  the  distribution  of  spectral  amplitudes  from 
test  to  test.  No  clear  patterns  in  this  data  are  evident,  suggesting  that 
because  of  the  distance  of  the  accelerometer  from  the  meshing  position  of 
the  faulty  gear,  the  fault's  effect  on  the  vibration  at  this  position  is 
attenuated.  It  might  be  noted  that  there  is  a  strong  spectral  component 
in  the  data  at  approximately  2600  Hz  and  at  harmonics  of  this  frequency. 

This  frequency  appears  to  be  the  frequency  of  passing  of  generator  commuta¬ 
tor  bars. 

Figures  54  and  55  show  spectra  of  acceleration  monitored  on  the  gearbox 
housing  for  the  case  of  good  and  bad  gears.  The  distribution  of  spectral 
amplitudes  shown  in  Figure  56  again  indicates  no  strong  correlation  of  the 
spectral  components  with  fault  presence.  This  might  again  be  expected  be¬ 
cause  of  the  number  of  interfaces  interposed  between  the  mesh  of  the  faulty 
gear  and  this  accelerometer  position. 

Mention  should  be  made  of  the  fact  that  in  all  but  one  of  the  tests  per¬ 
formed,  the  tooth  fault  was  present  on  only  one  tooth.  In  actual  gear 
operation,  it  is  quite  likely  that  the  type  of  fault  introduced  (pitch-line) 
pitting)  will  occur  on  a  number  of  teeth  and  in  most  cases  on  all  teeth. 

This  would  greatly  increase  the  forces  being  generated  in  the  gear  mesh 
which  would  likely  cause  stronger  manifestations  of  fault  presence  in 
acceleration  signals  monitored  at  positions  such  as  the  generator  output 
bearing  or  gearbox  housing.  In  this  study,  the  effects  of  faults  were 
strong  enough  at  the  input  bearing  position  to  provide  the  necessary  infor¬ 
mation. 
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FREQUENCY  (kHz) 

Figure  45.  Acceleration  Spectra  at  Input  Bearing  for  70$  Load. 
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O  GOOD  GEARS  SIDE  BAND  AMPLITUDES  IN  TERMS  OF 

A  GEARS  WITH  FAULTS  dB  DOWN  FROM  FIRST  HARMONIC 

AMPLITUDE 

Figure  47.  Distribution  of  Spectral  Amplitudes  at  Pinion  First  Harmoic  and  Sidebands 
at  Input  Bearing  Sensor  Location  Under  70$  Load. 
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Figure  48.  Distribution  of  Spectral  Amplitudes  at  Pinion  Second  Through  Seventh 
Harmonic  at  Input  Bearing  Sensor  Location  Under  70%  Load. 
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Figure  U9.  Acceleration  Spectra  at  Generator  Output  Bearing 
(Good  Gear) . 
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FREQUENCY (kHz) 

Figure  50.  Acceleration  Spectra  at  Generator  Output  Bearing 
(Gear  Fault). 
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Figure  51-  Distribution  of  Spectral  Amplitudes  at  Pinion  Meshing  Frequency  and  Sidebands 
at  Generator  Output  Bearing  Sensor  Location. 
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Figure  53.  Distribution  of  Spectral  Amplitudes  at  Pinion  Second  Through  Seventh 
Harmonic  at  Generator  Output  Bearing  Sensor  Location# 
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Tests  were  also  performed  to  determine  the  eflecto  of  speed  variations 
between  3700  and  U800  rpm.  While  there  was  some  variation  in  amplitudes  of 
the  meshing  frequency  and  its  harmonics  as  speed  was  varied,  no  large 
changes  in  the  amplitudes  as  the  result  of  torsional  resonances  were  ob¬ 
served. 

I 

STRAIN  GAGE  DATA 

— 

f 

Strain  gages  were  mounted  on  the  pinion  to  allow  monitoring  of  dynamic 
strain  while  the  gear  was  in  operation.  Figure  57  shows  typical  signals 
taken  from  a  gear  in  good  condition.  The  traces  represent  three  consecutive 
teeth  on  the  gear.  In  Figure  57a,  a  complete  meshing  cycle  may  be  observed. 
The  time  between  the  appearance  of  the  large  spikes  corresponds  to  the 
rotational  period  of  the  gear.  The  smaller  spikes  occur  at  the  gear  meshing 
frequency  and  represent  dynamic  hoop  stresses  of  the  gear  wheel.  The  higher 
frequency  component  evident  in  the  signal  is  caused  by  a  deficiency  of  the 
tape  recorder  which  allowed  some  of  the  carrier  frequency  to  "leak"  into 
the  data.  Figure  57b  presents  the  individual  tooth  contact  in  more  detail. 
Evident  in  this  picture  by  the  amplitude  of  strain  is  ;he  shift  between  one 
tooth-pair  and  two  tooth-pair  load  sharing.  As  may  also  be  seen  from  the 
figure,  shifts  in  load  sharing  are  not  as  abrupt  as  static  analysis  pre¬ 
dicts.  Based  on  the  figure,  a  contact  ratio  of  between  1.9  and  2.0  could 
be  found,  as  opposed  to  the  value  of  1.75  calculated  analytically.  This 
is  probably  the  result  of  dynamic  deflections  encountered  in  the  mesh. 

Figure  58  shows  torque  and  strain  data  from  a  gear  which  contained  a  small 
fault.  Figure  59  shows  similar  data  for  the  gear  containing  a  large  fault. 
Again  larger  strain  amplitudes  are  evident  at  the  gear  rotational  period 
while  smaller  oscillations  appear  at  the  gear-meshing  frequency.  There  are 
also  some  smaller  spikes  caused  by  interactions  of  the  dummy  gages.  The 
data  show  the  tooth  strains  to  be  significantly  different  for  faulty  gears. 
Oscillations  in  the  faulty  gears  occur  as  the  result  of  the  tooth-bending 
response  to  the  impulsive  loading  caused  by  the  fault.  There  is  also  some 
evidence  of  loss  of  contact  as  the  fault  passes  through  the  meshing  posi¬ 
tion. 


SUMMARY 


It  was  seen  that  the  effects  of  the  fault  introduced  into  the  geared  system 
may  be  observed  in  the  gearbox  vibration  time-domain  signals.  The  presence 
of  the  spikes  in  the  acceleration  signals  lends  credence  to  the  introduction 
of  an  impact  force  into  the  geared  system  as  a  result  of  the  fault.  It 
might  be  noted  here  that  impacts  in  the  mesh  might  also  be  caused  by  gear 
tooth  pitch  errors,  assembly  difficulties,  and  misalignment. 

Plots  of  the  input  torque  show  oscillations  about  the  mean  value.  This 
information  is  used  in  the  simulation  as  an  input  to  the  geared  system. 
Strain  gage  data  show  both  low  mesh  damping  and  indications  of  loss  of  con¬ 
tact  as  the  fault  passes  through  the  mesh.  Both  results  are  useful  as 
inputs  to  the  simulation  study. 
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(a) 


(b) 


Figure  57.  Strain  Fesponse  for  Good  Gears. 
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Strain  Signals  for  Gear  With  a  Small 
Fault . 
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Figure  59*  Strain  Signals  for  Gear  With  a 
Large  Fault. 
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The  primary  information  with  respect  to  discrimination  of  data  on  the 
basis  of  fault  presence  occurs  in  the  spectra  of  the  input  bearing  accel¬ 
eration.  Amplitudes  at  the  meshing  frequency  are  significantly  higher 
for  gears  with  faults,  and  a  sawtooth  consisting  of  multiples  of  shaft 
rotational  frequency  is  also  seen  for  faulty  gearing.  The  effects  of 
faults  on  sidebands  and  harmonics  are  not  readily  evident  since  distinc¬ 
tions  between  good  and  bad  gears  were  not  observed.  Response  information 
is  also  shown  to  be  a  strong  function  of  position  of  the  transducer  with 
respect  to  fault  position,  since  forces  generated  by  the  fault  tend  to  be 
attenuated  and  distorted  by  signal  transmission  paths. 
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SIMULATION  RESULTS 


The  equations  presented  in  the  Theoretical  Analysis  section  were  simulated 
and  solved  on  the  Ohio  State  IB4  System  370/165  computer.  Figure  60  pre¬ 
sents  a  flow  diagram  of  the  solution  process.  Stiffness  functions  for  both 
good  gears  and  those  with  faults  are  calculated,  and  Impact  amplitudes  due 
to  faults  are  computed.  These  data  are  then  used  in  the  torsional  dynamic 
simulation  where  the  geared  system's  time  response  is  calculated.  The  out¬ 
put  of  this  analysis  is  then  digitally  sampled  and  stored.  If  time-domain 
response  data  are  desired,  they  may  be  read  out  at  this  time.  The  data  may 
also  be  fed  through  a  Fast  Fourier  Transform  subroutine  which  calculates 
the  frequency  spectrum  of  the  system's  dynamic  response.  At  this  point, 
the  frequency  spectrum  may  be  either  read  out  directly  or  modified  by  any 
desired  transfer  function  and  then  read  out.  In  this  study,  output  data 
were  plotted  with  the  use  of  a  Calcomp  Plotter;  however,  the  dat*  are  also 
available  for  any  other  output  procedure  that  is  desired. 


MESH  STIFFNESS  FUNCTIONS 

As  a  first  step  in  the  analysis,  mesh  stiffness  functions  are  calculated. 
Figures  6l  and  62  present  the  calculation  procedure  for  good  gears.  In 
Figure  6l,  gear  compliance  is  plotted  versus  q  (distance  along  the  line  of 
action  from  the  pitch  point).  The  bending  compliance  for  the  gear  and  pin¬ 
ion  and  the  Hertzian  compliance  for  the  pair  of  teeth  are  shown  along  with 
the  total  tooth  pair  compliance  function.  The  use  of  load-sharing  consid¬ 
erations  allows  calculation  of  the  mesh  stiffness  function  apportioned  to 
each  tooth  on  the  gear  in  the  manner  shown  in  Figure  62.  The  parameter  KAP 
is  a  computational  dummy  variable  which  is  proportional  to  the  gear  rota¬ 
tional  position. 

Figures  63  and  64  represent  compliance  and  stiffness  functions  for  the  gear 
mesh  when  a  small  fault  (0.020  in.  wide)  is  present  on  a  tooth.  As  ex¬ 
plained  previously,  when  the  contact  strip  of  the  tooth  pair  is  contained 
within  the  fault  width,  contact  is  assumed  to  be  lost,  causing  Hertzian  com¬ 
pliance  and,  thus,  total  compliance  to  become  infinite.  This  makes  stiff¬ 
ness  values  zero  in  that  range,  as  shown  in  Figure  64.  In  the  torsional 
simulation,  mesh  damping  is  set  to  zero  when  the  stiffness  value  is  zero. 
Figures  65  and  66  show  similar  compliance  and  stiffness  curves  for  the  case 
of  a  fault  width  of  0.050  inch. 

One  other  stiffness  function  was  calculated  to  determine  the  effect  of 
change  in  load  sharing  a  system  response.  The  plot  of  Figure  67  shows  the 
stiffness  function  when  the  gear  mesh  contact  ratio  is  1.92. 


IMPACT  AMPLITUDE 


In  order  to  calculate  the  amplitude  of  the  torque  impulse  caused  by  a  tooth 
fault,  velocity  differences  as  a  function  of  fault  width  were  determined. 
Torsional  dynamic  simulation  was  used  for  this  purpose.  At  the  point  of 
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Figure  60.  Flow  Diagram  of  System  Simulation. 
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Figure  6l.  Tooth  Pair  Compliance  (Good  Gears). 


Figure  63.  Tooth  Pair  Compliance  (Small  Fault,  0.020 Inch  Wide 
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Figure  64,  Gear  Mesh  Stiffneus  "crsus  KAP  (Small  Fault, 
0.020  Inch  Wide) 
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Figure  65.  Tooth  Pair  Compliance  (Large  Fault,  0.050  Inch  Wide). 


initiation  of  contact,  the  fault,  mesh  stiffness,  and  damping  were  set 
to  zero,  allowing  the  geared  system  to  "freewheel",  i.e. ,  operate  with  a 
loss  of  contact  between  two  of  the  gears.  From  the  response  data  under 
this  condition,  Figure  68  shows  approximate  values  of  velocity  differences 
of  the  gear  and  pinion  and  time  spent  in  the  fault  as  a  function  of  fault 
width  on  the  tooth.  Fault  widths  are  scaled  in  terms  of  the  parameter  KAP. 
Based  on  the  desired  fault  width,  valves  may  be  read  from  the  plot  and  in¬ 
serted  into  the  torsional  dynamic  simulation. 


TIME-DOMAIN  RESFONSE 


Time-domain  response  data  was  one  of  the  outputs  available  from  the  simula¬ 
tion.  Figures  69-73  show  response  plots  of  ijq  ,  \|r?,  \|r3,  \Jr4,  and  for  the 
good  gears.  Amplitudes  are  presented  in  terms  of  radians  for  angular  dis¬ 
placement  and  rad/s2  for  angular  acceleration.  Evident  in  all  of  the  plots 
are  oscillations  at  the  gear-meshing  frequency  (the  result  of  the  mesh 
stiffness  variations)  and  also  at  the  input  shaft  rotational  frequency  (the 
result  of  the  input  torque  variation).  Figures  69  and  72  (tT  and  \|r4,  re¬ 
spectively)  show  less  of  the  effect  of  the  gear-meshing  frequency.  This  is 
to  be  expected  since  these  parameters  represent  relative  rotation  along 
shafting  of  the  system  rather  than  between  gears.  The  displacement  plots 
show  an  oscillation  about  a  nonzero  mean  level.  This  level  corresponds  to 
a  wrap-up  angle  of  the  torsional  system  as  the  result  of  transmitted  load 
torques.  The  acceleration  plot  of  Figure  73  shows  a  "beating"  at  the  shaft 
frequency  which  should  be  indicative  of  sidebanding  in  the  frequency  anal¬ 
ysis. 

Figures  74  and  75  represent  the  response  of  -y,  and  V?  when  the  large  fault 
(0.050  in.)  is  introduced  into  the  analysis.  The  fault  produces  the  large 
spike  in  the  response  with  a  repetition  rate  of  once  per  shaft  revolution 
which  corresponds  to  the  rate  at  which  the  fault  passes  through  the  meshing 
position.  The  relative  amplitude  of  the  spike  with  respect  to  the  steady- 
state  oscillation  is  larger  than  was  seer,  in  experimental  data.  This  is 
reasonable  since  there  is  likely  to  be  a  good  deal  of  attenuation  of  the 
spike  as  it  passes  from  the  gear  mesh  to  accelerometer  monitoring  positions 
on  the  gearbox. 


FREQUENCY  SPECTRA 


The  primary  comparison  of  experimental  and  theoretical  data  was  per¬ 
formed  in  the  frequency  domain.  The  frequency  transform  used  on  the  simu¬ 
lation  data  is  a  1024  input  point  FFT  analysis  performed  on  data  sampled 
at  20  kHz.  This  provides  output  information  fr<m  0-10  kHz  having  a  band¬ 
width  of  19.5  Hz,  which  is  comparable  to  that  obtained  with  the  real-time 
analyzer  on  the  experimental  data.  Six  averages  of  the  spectral  amplitudes 
are  performed  to  lower  noise  levels.  A  second-order  anti-aliasing  filter 
is  used  with  a  break  frequency  of  9  kHz*  Data  are  presented  from  0-7  kHz 
to  minimize  any  other  aliasing  effects  which  may  be  present. 
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Figures  76-80  show  frequency  spectra  for  1fit  i 'a#  i|r3,  ^4,  and  when 
no  gear  faults  are  present.  Obvious  peaks  in  the  spectra  may  be  seen  at 
3075  and  6150  Hz,  which  are  the  gear  meshing  frequency  and  its  first  har¬ 
monic,  respectively.  Also  present  are  sideband  frequencies  at  ±  75  Hz,  ± 
150  Hz,  etc.,  about  these  two  peaks.  A  large  peak  is  also  seen  at  75  Hz, 
the  input  shaft  rotational  frequency.  There  are  also  frequency  components 
in  the  regions  near  1500  and  4500  Hz.  The  exact  origin  of  these  components 
has  not  been  ascertained,  but  they  appear  in  the  experimental  results  and 
may  be  the  result  of  nonlinear  effects  of  the  system.  These  peaks  are  more 
predominant  in  the  acceleration  data  of  Figure  80.  Torsional  resonant  fre¬ 
quencies  for  the  simulation  were  calculated  using  a  Holzer  analysis,  the  re 
suits  of  which  are  shown  in  Table  IV.  Since  the  mesh  stiffness  is  variable 


TABLE  IV.  SIMULATION  TORSIONAL  RESONANCES 
VIA  HOLZER  ANALYSIS 


One  Tooth- Bair 

Two  Tooth-Pair 

Stiffness 

Stiffness 

(Hz) 

(Hz) 

458 

464 

918 

943 

4247 

5644 

5264 

6974 

mean  values  of  single  tooth-pair  stiffness  and  double  tooth-pair  stiffness 
were  used  in  calculating  two  sets  of  natural  frequencies.  Actual  resonant 
values  will  lie  somewhere  within  this  range;  however,  no  strong  components 
in  the  spectra  were  observed  at  these  frequencies. 

One  fact  which  may  be  noted  from  the  spectra  is  that  the  relationship  of 


h  =  u>2(\jr?)  (57) 

is  not  followed  at  all  frequencies.  It  is  felt  that  thi-  is  a  result  of 
the  smoothing  and  computational  errors  induced  >y  the  double  integration 
used  to  calculate  angular  displacement  from  acceleration.  However,  the 
amplitudes  of  the  frequencies  at  which  differences  occur  are  such  that 
lictle  effect  on  the  meshing  frequency  is  expected. 

Figures  8l-84  show  the  frequency  spectra  of  and  \|r2  when  fault  effects 
have  been  introduced  into  the  analysis ;  \jf2  and  if2  have  been  analyzed  since 
they  represent  the  relative  displacement  and  acceleration  of  the  gear  mesh 
in  which  the  fault  is  encountered  and,  thus,  should  show  the  strongest 
effects  of  the  faults.  In  Figures  8l  and  82,  the  0.020 -inch -wide  fault  has 
been  introduced,  while  in  Figures  83  and  84  the  0,050-inch -wide  fault  was 
introduced.  In  the  case  where  faults  are  introduced,  amplitudes  at  the 
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meshing  frequency  and  its  first  harmonic  increase  as  shown  in  Table  V, 


The 


TABLE  V.  SPECTRAL  AMPLITUDES  OF  SIMULATION  RESULTS  (dB) 


Frequencies 

Good 

Gear 

Small 

Fault 

Large 

Fault 

Good 

Gear 

Small 

Fault 

Large 

Fault 

Gear  Meshing 

31.3 

35.5 

37.1 

42.5 

46.7 

48.4 

1st  Lower  S.B. 

25.3 

29.5 

31.4 

34.4 

40.7 

42.4 

2nd  Lower  S.B. 

8.4 

15.3 

19.4 

7.8 

23.0 

30.0 

1st  Upper  S.B. 

27.0 

23.8 

23.0 

41.3 

36.4 

35.0 

2nd  Upper  S.B. 

1 6.6 

12.7 

17.2 

25.8 

24.1 

29.2 

1st  Harmonic 

28.7 

33.3 

35.2 

57.2 

57.9 

58.6 

1st  Lower  S.B. 

27.7 

28.2 

27.8 

52.3 

50.3 

51.1 

2nd  Lower  S.B. 

20.1 

18.0 

22.7 

37.5 

42.9 

45.7 

1st  Upper  S.B. 

25.6 

27.4 

23.3 

48.1 

48.1 

47.0 

2nd  Upper  S.B. 

17.1 

20.1 

16.7 

43.1 

39.2 

40.5 

increases  are  greater  at  the  meshing  frequency,  and  a  further  increase  is 
seen  with  an  increase  in  fault  size.  The  changes  in  sideband  amplitudes  do 
not  follow  a  distinct  pattern.  The  changes  in  amplitude  at  the  meshing  fre¬ 
quency  are  not  as  great  as  those  seen  in  the  test  data,  but  the  simulation 
response  data  is  quite  repeatable,  at  least  with  respect  to  peak  amplitudes, 
so  the  statistical  variations  encountered  in  the  experimental  data  are  not 
present.  Comparison  of  the  figures  also  shows  the  introduction  of  the  saw¬ 
tooth  effect  of  peaks  at  multiples  of  the  shaft  rotational  frequency.  Ampli¬ 
tudes  of  this  phenomenon  are  further  increased  when  fault  size  is  increased. 
An  interesting  pattern  may  be  observed  in  these  peaks.  At  intervals  of 
approximately  450  Hz,  a  frequency  beating  occurs  in  which  the  envelope  of 
amplitudes  alternately  increases  and  decreases.  This  spacing  corresponds 
to  the  lowest  torsional  resonance  of  the  system.  It  is  interesting  to  note 
that  a  similar  effect  was  seen  in  the  experimental  data  illustrated  in 
Figure  39, 

An  experimental  transfer  function  between  acceleration  at  the  gear  mesh  and 
gearbox  input  bearing  acceleration  was  obtained  as  described  in  Appendix  V, 
This  information  was  then  used  along  with  simulation  response  data  to  pro¬ 
duce  the  frequency  spectra  of  Figures  85,  86,  and  87,  which  represent  accel¬ 
eration  of  a  good  gear,  a  gear  with  a  0. 020-inch-wide  fault,  and  a  gear  with 
a  0.050-inch-wide  fault,  respectively.  While  amplitude  increases  at  the 
meshing  frequency  and  the  sawtooth  effect  is  somewhat  apparent,  overall  cor¬ 
relation  with  experimental  results  is  not  extremely  good.  Peak  spectral 
amplitudes  of  the  acceleration  at  gear  meshing  frequency  and  its  first  har¬ 
monic  fall  in  the  range  of  5-10g,  which  is  the  same  range  noted  in  experi¬ 
mental  data.  However,  differences  are  observed  in  the  overall  shapes  of  the 
spectra.  One  reason  for  the  differences  is  the  lack  of  inclusion  of  bending 
response  in  the  simulation,  since  it  is  purely  torsional  in  nature.  The 
manner  of  producing  the  transfer  function  is  also  suspect.  Much  more  work 
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is  necessary  in  order  to  develop  techniques  to  obtain  more  reliable 
transfer  function  information  of  gear  mesh  systems  of  the  type  included 
in  this  study. 

Information  was  desired  on  the  effect  of  changing  certain  system  input 
parameters  on  the  frequency  spectra.  Figures  88  and  89  show  the  effect 
of  changing  the  amount  of  torque  oscillation  while  keeping  the  mean  values 
constant.  In  Figure  88  there  is  no  oscillation,  and  in  Figure  89  the  oscil¬ 
lation  is  approximately  23 io  of  the  mean  value.  This  contrasts  to  the 
oscillation  of  approximately  4 6$  which  was  used  in  the  previous  results 
and  is  shown  in  Figure  77.  As  seen  from  these  figures  there  is  a  change 
of  one  or  two  decibels  in  the  amplitudes  at  the  meshing  frequency  and  its 
first  haimnic.  A  dramatic  change  is  evident,  though  in  the  sideband 
frequencies.  The  amplitude  and  number  of  sidebands  decrease  as  a  function 
of  the  torque  oscillation.  This  is  in  opposition  to  the  lack  of  any 
dramatic  changes  in  sidebanding  as  a  result  of  fault  introduction  in  either 
experimental  or  simulation  spectra.  This  would  lead  to  the  conclusion  that 
fault  presence  is  not  the  major  mechanism  involved  in  changes  in  sideband¬ 
ing. 

The  two  mechanisms  introduced  into  the  simulation  as  a  result  of  fault 
presence  were  a  special  stiffness  function  for  a  faulty  tooth  and  a  torque 
impact  due  to  the  fault.  The  relative  effects  of  these  two  mechanisms  are 
shown  in  Figur  -  90  and  91.  In  Figure  90,  an  impact  is  introduced  without 
the  special  stiffness  function;  the  reverse  is  true  in  Figure  91-  Both 
spectra  show  similar  characteristics  of  relatively  the  same  order  of  mag¬ 
nitude  and,  thus,  are  responsible  in  combination  for  producing  the  overall 
response  to  a  faulty  gear  as  seen  in  Figure  75* 

Figure  92  shows  a  spectrum  for  the  system  for  good  gears  when  the  contact 
ratio  of  1.92  was  used.  Amplitudes  are  slightly  lower,  as  might  be  expected, 
since  the  overall  effect  is  to  make  the  system  somewhat  stiffer  and  because 
the  length  of  time  spent  in  single  tooth-pair  contact  is  shorter. 

As  mentioned  in  the  section  on  experimental  results,  it  is  unlikely  that  a 
tooth  fault  of  the  type  used  in  this  study  would  occur  on  only  one  tooth. 

For  this  reason,  the  spectrum  of  Figure  93  was  run  showing  the  effect  of 
the  fault  stiffness  variation  when  a  fault  was  introduced  on  every  tooth. 

As  may  be  seen,  amplitudes  at  all  frequencies  are  greatly  increased.  There 
is  no  evidence  of  the  sawtooth  effect  since  the  excitations  due  to  the 
fault  occur  at  the  gear-meshing  frequency.  The  effect  of  the  torque  impact 
is  not  shown  since  an  impact  occurring  at  a  repetition  rate  corresponding 
to  the  gear-meshing  frequency  would  produce  significant  energy  at  fre¬ 
quencies  beyond  the  present  capabilities  of  the  model.  Modifications  could 
produce  this  capability  but  at  extreme  penalties  in  computational  costs. 


SUMMARY 


Gome  evidence  of  the  fault  in  the  form  of  spikes  was  noted  in  the  time- 
domain  response.  Changes  at  the  gear-mesh  frequency  and  harmonics  are  also 


127 


or  09 —  oroh  or  02  oro  ao-oz- 

ibq)  aflnintiwkMwn 


i  < 

-or  Oh- 


128 


I  I  I  I  I 

00'  09  00"  Oh  Off"  00  00"  0  DO"  02" 

(80J  30/1 II  “WHIM)  01 


I  o 

oroh- 


s> 

*H 


130 


Lde  Fai 


132 


133 


FREQ 


_ _ 

6 


evident  in  the  frequency  spectra.  When  a  fault  is  introduced,  the  ampli¬ 
tude  of  the  gear-meshing  frequency  increases,  a  saw tooth  effect  at  multiples 
of  the  input  shaft  frequency  occurs,  tnd  a  periodicity  in  the  sawtooth 
corresponding  to  approximately  the  first  torsional  resonance  of  the  dynamic 
model  is  seen.  All  three  effects  correspond  to  changes  noted  in  experimen¬ 
tal  data. 

It  was  found  that  changes  in  sidebanaing  occurred  as  a  result  of  an  input 
torque  oscillation  rather  than  as  a  direct  result  of  fault  introduction. 

It  was  also  determined  that  changes  in  frequency  spectra  of  gears  with 
faults  were  a  composite  effect  of  stiffness  variations  and  torque  impacts 
caused  by  the  faults.  This  would  lead  to  the  assumption  that  similar  spec¬ 
tral  changes  would  be  noted  when  types  of  faults  which  do  not  cause  loss 
of  contact  and  subsequent  stiffness  variations  are  analyzed.  Such  faults 
(e.g. ,  wearing  or  scoring)  could  induce  impacts  alone  into  the  geared 
system.  Finally,  when  faults  are  present  on  each  tooth  of  the  gear,  the 
amplitudes  of  both  time-domain  and  spectral  information  were  significantly 
increased. 


CONCLUSIONS 


As  a  general  conclusion  of  the  research  program,  it  may  be  stated  that  the 
analysis  and  modelling  procedures  employed  provide  an  effective  method  for 
determining  the  changes  encountered  in  a  geared  system's  response  as  a  re¬ 
sult  of  tooth  surface  faults.  In  analyzing  the  comparison  of  response  data 
for  good  and  bad  gears  from  both  the  simulation  and  experiments,  the  follow¬ 
ing  more  specific  conclusions  are  offered: 

1.  Fault  presence  in  a  single  tooth  of  a  gear  causes  three  major 
differences  in  system  frequency  spectra;  i.e., 

a.  increases  in  amplitude  at  the  gear-meshing  frequency, 

b.  a  sawtooth  effect  consisting  of  energy  at  multiples  of 
the  input  shaft  frequency,  and 

c.  a  periodicity  in  the  sawtooth  effect  corresponding  to  the 
system's  lowest  torsional  resonance. 

2.  The  effects  of  fault  presence  on  harmonics  of  the  gear-meshing 
frequency  and  sidebands  of  both  the  gear-meshing  frequency  and  its 
harmonics  do  not  follow  any  distinct  pattern. 

3.  Fault  presence  may  be  noted  in  time-domain  data  by  the  presence  of 
spikes  in  the  data. 

Simulation  response  data  provided  the  following  further  conclusions: 

1.  The  major  factor  affecting  sideband  generation,  in  this  study,  was 
an  input  torque  oscillation  which  was  felt  to  be  the  result  of 
mounting  fixtures  and  system  alignment. 

2.  Fault  effects  were  the  results  of  both  stiffness  variations  (the 
result  of  loss  of  contact  at  the  fault)  and  torque  impacts  (the 
result  of  velocity  differences  when  contact  is  reestablished) . 

3.  The  relative  effects  due  to  the  stiffness  variations  and  torque 
impacts  are  of  the  same  magnitude. 

A  deficiency  in  the  model  was  evidenced  by  difficulties  encountered  in  pre¬ 
dicting  the  exact  amplitudes  (in  either  the  frequency  or  time-domains)  of 
vibrations  as  measured  on  the  outside  of  the  gearbox.  Reasons  for  this 
were  as  follows: 

1.  The  model  analyzed  purely  torsional  motion  without  accounting  for 
bending  effects. 

2.  Accurate  transfer  function  information  between  the  gear  mesh  and 
accelerometer  positions  was  not  obtained. 

While  the  study  did  not  develop  a  specific  vibration  diagnostic  technique, 
it  is  felt  that  the  results  obtained  will  be  useful  in  the  following  ways: 

1.  An  understanding  and  explanation  (which  were  not  previously 

available)  of  the  mechanisms  involved  in  the  gear  mesh  as  a  result 
of  fault  presence  have  been  developed. 
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2.  A  simulation  of  a  geared  system's  response  to  one  type  of  tooth 
surface  fault  has  been  developed,  spectral  changes  as  a  result  of 
fault  presence  have  been  noted,  and  the  results  have  correlated 
with  experimental  results. 

3.  Ghe  model,  as  developed,  shows  distinct  possibilities  for  exten¬ 
sion  to  other  fault  types  to  allow  its  use  as  a  general-purpose 
gear  response  predictor. 
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RECOMMENDATIONS 


Based  on  the  results  of  this  study,  there  are  two  areas  in  which  further 
development  work  should  be  considered.  One  involves  extension  of  the 
capabilities  of  the  model  as  developed,  and  the  other  involves  more  sophis¬ 
ticated  analytical  work  which  may  be  performed  in  the  area  of  gear  dynamic 
analysis. 

Further  development  of  the  model  would  consist  of: 

1.  Testing  to  provide  a  larger  data  base  for  analysis  of  the 
model,  including 

a.  more  gears, 

b.  more  types  of  faults, 

c.  variation  in  fault  position,  and 

d.  variations  in  the  severity  of  the  faults. 

2.  Testing  to  provide  better  transfer  function  information 
between  the  gear  mesh  and  external  gearbox  transducer 
positions. 

3.  Combination  of  bending  effects  with  simulation  results. 

4.  Redesign  of  test  stand  fixtures  to  eliminate  the  torque 
oscillations  (if  possible). 

5.  Addition  of  more  lumps  to  the  torsional  model;  i.e., 
providing  more  detailed  analysis  of  the  driving  section 
of  the  test  stand  and  of  the  generator. 

The  acquisition  of  more  data  could  provide  better  statistical  information 
on  data  variations.  The  amount  of  data  taken  in  the  study  was  not  large 
enought  for  a  very  valid  statistical  analysis.  The  use  of  better  transfer 
function  information  and  the  inclusion  of  bending  effects  should  permit 
better  correlation  of  specific  amplitudes  between  the  simulation  results 
and  measured  experimental  data.  The  torque  variation  was  the  major  factor 
in  changing  sideband  amplitudes.  If  this  variation  is  diminished  or  elimi¬ 
nated,  more  subtle  changes  as  the  result  of  fault  presence  may  be  distin¬ 
guishable.  The  addition  of  more  lumps  to  the  torsional  model  is  a  natural 
extension  which  brings  the  model  closer  to  its  natural  continuous  state. 
Limitations  axe  reached,  however,  at  the  point  where  increased  computational 
burdens  are  not  offset  by  improvements  in  the  simulation  response. 

There  axe  two  axeas  in  which  further  theoretical  development  could  proceed; 

i.e., 

1.  The  inclusion  of  provisions  to  allow  tooth  contact  off 
the  line-of-action. 

2.  Two-dimensional  modelling  of  the  geax  tooth. 
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If  analysis  were  extended  to  motion  off  the  line-of-action,  noninvolute 
contact  could  be  considered.  This  would  provide  results  for  faults  which 
cause  major  tooth  profile  modifications  such  as  wear  or  manufacturing 
errors.  It  could  also  allow  the  inclusion  in  the  analysis  of  noninvolute 
gears.  The  inclusion  of  another  dimension  in  the  analysis  would  allow 
results  to  apply  to  spur  gears  with  larger  facewidths,  as  well  as  bevel 
and  helical  gears.  There  would  also  be  provisions  for  analyzing  faults 
which  are  more  two-dimensional  in  nature  them  pitch-line  pitting;  e.g. , 
scoring,  spalling,  etc. 

As  a  general  recommendation  for  both  extension  of  the  present  model  and 
the  development  of  more  sophisticated  models,  hybrid  computation  should 
be  considered.  The  complexity  of  the  model,  switching  facilities,  logic, 
and  function  storage  required  make  pure  analog  simulation  unfeasible,  but 
purely  digital  analysis  also  presents  problems.  The  large  frequency  range 
which  the  system  exhibits  requires  very  small  computing  intervals  for  any 
numerical  integration  scheme,  making  computing  quite  expensive.  Thus,  a 
hybrid  simulation  which  would  allow  analog  solution  of  the  differential 
equations  with  the  digital  capabilities  of  logic,  switching,  and  iterative 
computations  appears  to  be  the  optimum  computation  technique. 
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APPENDIX  I 
GEAR  MEASUREMENTS 


Typical  measurements  for  the  gears  used  in  the  test  program  were  obtained 
from  the  manufacturer  and  are  shown  in  Figures  9h,  95,  and  96.  Included 
in  the  charts  are  profile  and  spacing  checks  as  well  as  a  lead  check  on  the 
pinion.  From  Figures  9*+  and  96,  it  may  be  seen  that  the  spur  pinion  and 
generator  drive  gear  are  very  close  to  a  perfect  involute  with  variations 
of  less  than  0.0001  inch.  The  chart  of  the  idler  gear  (Figure  95)  shows 
the  presence  of  profile  modifications  of  up  to  0.0002  inch.  These  modifi¬ 
cations  are  of  a  type  used  to  compensate  for  gear  tooth  deflections  under 
load.  The  spacing  charts  for  the  gears  show  a  maximum  tooth-to-tooth 
variation  of  0.0001  inch  and  total  spacing  error  of  0,0005  inch.  No 
significant  lead  error  on  the  pinion  is  indicated  by  the  lead  check. 

Figures  97,  98,  and  99  show  profile  checks  measured  after  testing  was  com¬ 
pleted  for  several  of  the  gears.  There  is  little  indication  of  any  large- 
scale  wear  on  the  gears,  which  would  be  the  expected  result  because  of  the 
limited  number  of  hours  of  service  of  the  gears.  The  charts  in  Figures  98 
and  99  show  the  gears  which  had  faults  manufactured  into  them.  The  depths 
of  the  faults  are  not  accurately  portrayed  because  of  limitations  of  the 
measuring  device.  There  is  no  indication  from  the  measurements  that  there 
were  any  significant  profile  variations  from  tooth  to  tooth  on  the  same 
gear  or  from  gear  to  gear.  There  is  some  wear  (up  to  0.0001  in.)  shown 
near  the  tooth  roots  on  some  of  the  gears,  but  this  is  not  too  significant 
and  falls  within  manufacturing  specifications. 
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Figure  96.  Gear  Measurements  for  Typical  New  Generator  Drive  Gear. 
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APPENDIX  IT 

COMFUTER  PROGRAM:  TORSIONAL  DYNAMIC  SIMULATION 


The  torsional  dynamic  equations  are  solved  via  the  IBM  Continuous  Systems 
Modelling  Program  (CSMP) .  The  response  data  are  sampled  and  stored,  and 
then  a  frequency  analysis  is  performed  with  the  use  of  the  F0UR2  (Fast 
Fourier  Transform)  Subroutine  of  the  IBM  Scientific  Subroutine  Package. 


The  following  input  data  are  required  on  PARAM  cards  for  the  dynamic  anal¬ 
ysis: 


GOOD 

WDTH 

PI 

NP 

J1 

J2 

J3 

J4 

J5 

B12 

B23 

B34 

B45 

K12 

K45 

R2 

R3 

r4 

SSHAFT 

KG 

KT 

RG 


indicator  of  whether  geared  system  has  any  faults.  If  GOOD  >  2, 
system  is  in  "good"  condition, 
fault  width,  rad 

3.14159 

number  of  teeth  on  the  pinion 

inertia  of  shafting  and  input  gear,  in.-lb*s 

inertia  of  spur  gear  pinion,  in.-lb*s2 

inertia  of  idler  spur  gear,  in.-lb*s2 

inertia  of  spur  gear,  in.-P  *s2 

inertia  of  generator,  in.-lb*s2 

torsional  damping  coefficient  of  input  shafting,  in.-lb/rad* s -l 

linear  damping  coefficient  of  pinion-idler  mesh,  lb/in. 's-1 

linear  damping  coefficient  of  idler-gear  mesh,  lb/in. ’s1 

torsional  damping  coefficient  of  spur  gear  shaft,  in. -lb/rad- s _1 

torsional  stiffness  of  input  shafting,  in.-lb/rad 

torsional  stiffness  of  spur  gear  shaft,  in.-lb/rad 

spinion  base  radius,  in. 

idler  base  radius,  in. 

gear  base  radius,  in. 

input  shaft  speed,  rpm 

generator  voltage  constant,  v/rad*s~l 

generator  current  constant,  in.-lb/A 

load  resistance,  fi 


To  increase  computing  efficiency,  a  set  of  parameters  is  input  which  starts 
computed  variables  at  a  point  in  time  at  which  the  system  is  in  steady- 
state  oscillation.  These  parameters  eliminate  the  computational  transients 
encountered  when  starting  ?n  analysis  from  zero. 


PS  ID 

= 

initial 

value 

of 

Vi  J 

rad 

PS2D 

initial 

value 

of 

'\fs  J 

rad 

PS3D 

= 

initial 

value 

of 

t3> 

rad 

ps4d 

initial 

value 

of 

^4  > 

rad 

re  iv 

initial 

value 

of 

fl  > 

rad/s 

PS2V 

initial 

value 

of 

rad/s 

PS3V 

=: 

initial 

value 

of 

^3  ) 

rad/s 

ps4v 

= 

initial 

value 

of 

^4J 

rad/s 

TH2D 

initial 

value 

of 

rad 

El 

initial 

value 

of 

derivative  of 

E2 

initial 

value 

of 

filter  output 

lr-0 


here  are  ten  parameters  which  are  required  by  sampling  and  storage  rou¬ 
tine  and  the  FFT  analysis: 


NPTS 

INCH 


IDATA 

AVG 

STIME 

WAA 

STIMEI 

ISIGN 

IFOM 

NDIM 


number  of  points  in  FFT  analysis  [must  be  power  of  2  and 
equal  to  TABLE  NR  (l)  'j 

Sampling  increment,  s  (less  than  or  equal  to  the  inverse  of 
twice  the  highest  frequency  of  interest.  It  must  also  be 
an  integer  multiple  of  TIMER  parameter  BELT) . 
number  of  Fourier  components  to  be  printed  out. 
number  of  averages  of  Fourier  components  to  be  performed, 
starting  time  of  sampling,  procedure,  s 
break  frequency  of  anti-aliasing  filter,  rad/s 
STIME 
-1 
1 
1 


Storage  declarations  are  required  based  on  the  number  of  points  in  the 
transform  (NPTS)  and  the  number  of  averages  to  be  performed  (AVG.).  The 
parameters  are  dimensioned  in  the  following  manner  with  actual  numbers  re¬ 
placing  the  symbols: 

W(NPTS),  AMP  (NPTS),  and  DBAMP  (NPTS). 

DATA  (NPTS),  RSPN  (NPTS),  and  FRC  (NPTS,  AVG). 

The  format  of  the  parameter  cards  is: 

Column  1-5 :  PARAM 
Colum  6 :  Blank 

Column  7-72:  Parameter  name  =  (any  format),  etc. 

Example:  PARAM  K12  =  8.43E  +  04,  K45  =  5.0E  +  05 

There  are  several  initial  condition  parameters  required  also  : 

TABLE  NR  (1)  =  NPTS 

INCON  N  =  0,  NI  =  1,  I  =  1,  FLT  =  0 

INCON  K23  =  initial  value  of  pinion-idler  mesh  stiffness 

INCON  K34  =  initial  value  of  idler-gear  mesh  stiffness. 

The  arbitrary  function  generator  option  (AFGEN)  is  used  to  input  the  gear 
mesh  stiffness  functions.  The  stiffness  calculation  programs  provide 
punched  card  output  formatted  as  required  for  the  torsional  dynamic  pro¬ 
gram. 

The  input  to  the  dynamic  (DYNAM)  portion  of  the  program  is  the  driving 
torque  (TQIN)  of  the  geared  system.  Any  constant  value  or  constant  value 
plus  oscillation  may  be  used  depending  upon  measured  results  from  the  test 
stand. 


The  parameters  of  the  program  corrr  r~md  to  the  parameters  of  the  theoreti¬ 
cal  analysis  as  follows: 
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> 

II 

|S> 

TH1D0T 

=  01 

A2  =  02 

TH2D0T 

• 

=  0 

TH2  =  02 

A3  =  Ra 
r2 

03  TH3D0T 

=  03 

a4  =  04 

thUdot 

=  04 

A5  =  05 

TH5D0T 

=  05 

PS1  =  ilq 

PS1D0T  = 

k 

PS1DDT  =  ilfi 

PS2  =  i|r2 

PS2D0T  = 

PS2DDT  =  i )r> 

K 

U> 

11 

03 

PS3D0T  = 

• 

PS3DDT  =  \|r3 

PS4  = 

PS4D0T  = 

t4 

PS4DDT  =  ilr4 

The  statement  : 

E2DT  =  (WAA**2)  *PS2-2.*.6*WAA*EDT-(WAA**2)*E 

represents  a  second-order  anti-aliasing  filter  with  break  frequency  of  WAA. 

The  NOSORT  section  of  the  program  picks  the  value  of  gear  mesh  stiffness 
needed,  determines  if  the  system  is  to  be  analyzed  with  a  fault,  and  picks 
the  special  stiffness  and  torque  impact  as  may  be  required.  The  gear 
stiffness  is  clocked  on  a  parameter  going  from  0  to  1  for  start  to  end  of 
the  tooth.  The  program  next  routes  the  computer  through  the  fault  modifi¬ 
cation  portion  of  the  program  depending  upon  the  value  of  GOOD. 

The  first  portion  of  the  fault  analysis  determines  if  a  faulty  tooth  is  in 
mesh  and,  if  so,  uses  the  modified  stiffness  function.  The  program  next 
sets  mesh  damping  to  zero  whenever  the  stiffness  is  zero.  The  torque  im¬ 
pact  is  also  introduced  at  the  appropriate  contact  point.  Three  parameters 
are  associated  with  the  fault  insertion: 

VPDEL  =  velocity  difference  of  pinion  caused  by  fault 
VGDEL  =  velocity  difference  of  gear  caused  by  fault 
FLTAU  =  time  of  impact  duration 

The  following  section  of  the  program  forms  the  sampler.  The  desired  data 
point  is  sampled  at  the  sampling  interval  (INCR)  and  stored  in  an  array  of 
dimension  NFTS  X  AVG. 

When  the  dynamic  analysis  is  complete,  the  program  enters  the  TERMINAL 
segment  to  perform  the  Fourier  Transform  and  output  the  uata.  The  routine 
performs  successive  transforms  on  the  data  and  then  uses  an  arithmetic 
averaging  of  the  coefficients.  The  amplitudes  of  the  transformed 
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coefficients  are  converted  to  decibels  referenced  to  whatever  value  is 
convenient.  The  frequency  response  data  are  then  available  on  a  set  of 
cards  containing  the  frequency  and  decibel  amplitudes.  Cards  are  also 
punched  with  time  and  time-domain  amplitude.  The  data  are  available  for 
whatever  further  processing  or  plotting  is  desired.  A  listing  of  the  pro¬ 
gram  follows,  with  a  sample  set  of  parameters  used  in  one  of  the  runs. 
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llTlfc  IDRSIUNAL  DYNAMIC  SIMULATION  TIE  HI  I  S.' I  uEAKbUX 

• 

STORAGE  M« ( 1  ) 

F  I  XFI)  N.NOIM,  IS  FAN,  I  FORM,  NR,  I.NP.NPTS,  I  OAT  A , M, M I , Y , 01 , N I , A VO, J 
/  DIMENSION  W<  10241,  AMP{  1024), D8AMPM024) 

/  COMPLEX  0ATAM024I.RSPN<  10?4), CARS. REAL, FRCl  1024,6) 

TARIF  NW C 1  I  —  1024 
I NCON  N=0, N I = 1 
1NC0N  823=7.52 

I  NCON  K23  =  ll.  t  i-'.E+05,K34=8.599E+0:> 

I NCON  1=1 
I NCON  FL T  =0. 

PARAM  0000=1 
INIT 

I  WAA=2.*P  l*FAA 

wi=?.#p!«sshaf  r/6o. 

T  H 1  V  =  W  1 

ATTH=2.«P|/A1. 

F L  1 C  A L  =  (  J2*VPI)EL+J38VGDEL  )  /  (  1  2.  *HT  All#  (  R2  +  R3  )» 

NOSORT 

DO  5  I=1«NPTS 
DATA! 11=0. 0 
5  CONTINUE 

SORT 

*  THF  AFOFN  DETERMINES  THE  OEAR  STIEENESS  AS  A  FUNCTION 

*  OF  OEAR  ROTATION.  STENS  IS  EUR  A  GOOD  OEAR.  STFNS?  IS 

*  FOR  A  GEAR  WITH  A  FAULT. 

AFOFN  STFNS  =  0.0  , 0. H 0  1  A, 0. 0373 . 0.  8652 , 0. 07 AS , 0. 8684, .  .  . 

0.1115,0.8709,0. 1A«5, 0.8  728,0.  1853,0.4  741,0.2221,0.8748.... 
0.  2587 , 0.  8749,  0.  2952,0.4  74  3, 0.  3316, 0. 8732,0. 3679,0.8  71  4,  ..  . 
0.4040,0. 8690,  0.4401, 0.4  72''.  0.4 760, 0. 4  7R2 , 0. 5 11 H , 0.4834 , . . . 
0. 5475,0.4879.0. 5831 ,0.4918,0.6185,0.4951 ,0.6538,0.4978, .. . 
0.6890,0. 794 «, 0.7241 ,0.80A6,0. 7890, 0.M1 37 ,0. 7938,0.8222, . . . 
0.828  5,0.8  30' ,0.  8631,0.8373,0.  8978,0.  44 39 , 0. 9 3 1 8 , 0. 8499  ,  .  . . 
0.9660,0.8552,  1.0000,0.8599 

AFOFN  ST  F'JS2  =  0.  0  ,  0.  8568, 0.  037  3, 0.  H603,  0.  0745, 0.  8631 ,..  . 

0. 1 1 1 8, 0. H65S, 0. 1 485, 0.86  7 l , 0. 1 85  3, 0. 8682,0. 2221 ,0.8688 ,.. . 
0.  258  7,0.8  645, 0.  2 ‘•'52, 0.  86  /  8 , 0.  3  3  16,  0.  466  3,0.  3  679 , 0 . 8  64  3  ,  .  .  . 
0.4  04  0,  0.  86  W,  0.4401 , 0.44  v*.  n.4 /M),  0.0000,0.5)  1  8,0.0000  .... 
0.  54  75, 0. 0000,  0.  5»'31 , 0.  0000,  0,  6]  85  ,  O.  0000, 0.  65  34  ,  O.  495  3  ,  .  .  . 
0.6890, 0.  7923, 0.7241 ,0.8018, 0. 7590, 0. 4108 . 0. 7934,0,8191  ,.. . 
0.8  285, 0. 8  268, O.  86 3 1 , 0. 8 X  3 8 , 0.  8975,0, 840 1 , 0. 9 3 1 8 , O.  8460,  . .  . 
0.9660,0. H5 1  1  ,  1 .0000,0.8  554 

* 

DYNAM 

8 

«  TORSIONAL  dynamics  OF  SYSTEM 

T0IN=330. +150. *S!N(wi*T!Mt ) 

A1  =  IK128PS1*H)2*PS1I'0T  +  H>|N»/J1 

A2  =  (K238(  R2**2  )-PS2i  R23M  R2<=*2  )  <‘PS?DOT-K  12«PS  1-81  2*PS1  DOT..  . 
-ELT8R2 I / J  2 

A  3  —  ( K348R3-P4-PS  3  +  R348R  38R48PS 3UUT-K2  38R38K28PS2. .  . 

-B23=>H  38R28PS2IMI1  +  EL1*R?  l»R3/(  J3*R2) 

A4  =  l  K45*PS4  +  R4  58PS4DUT-K34  >  (  R 4*<=  2  )  <'  PS 3-834*  (  R48*2  )  *PS3I)0T  )  /.|4 

A5  =  -(K45*PS4  +  R4  5*PS4()UT+K  T*  IGT/J5 

I0  =  K081H50(II/PG 

PS  1  DOT  =  A  2- A l 

PS2D01 = A  3- A  2 

PS3DDT  = A4-A  3 

PS400T  =  A5- A4 

PSinnT  =  iNK;»L i rsi v.psiddi  i 
PS2DUT  = I NTORL ( PS2V, PS20D1 ) 
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PS1l)!ll-|Mli.“L  !  PS  AV  .  PS  3001  ) 

PS400I  -  I  MTt'.w  r  I  P  '>  4  V «  P  S4DDI  > 

PS1  =  !NT<;i<l  (  PS  1 1', MS  1  HlJ  1  ) 
ps?=  in  1  gri  ( pspu. pspoih  j 
psi= ini  gri.  i  psao.ps  iooti 

PS4=|N1(;*>L(  PS4I),  PS4IX1T) 
iHionTsiNir.miMi.Ai ) 

th2dot=ps  i  ut  it  ♦  th  i  in  ii 

TH300T  =  I  PS2U01  *  IH^iuiT  )*R2/R3 
T  H4linl  =PS3DOT+R 333^300 T/w4 
TH50(IT  =  PS40nT  +  lH4niiT 
TH?=  INTGP  L  (  TH?|),  TH?U(II  | 

*  ANT  I  -AL  !  AS  I  MG  F  I  L  T  (- « 

F?OT=-2.  *.  6*WAA*FI)l-{  WAA**2  )*F+(  WAA**2  )*PS2 
POT  = INTGRL  m  ,F?Dl  ) 

MINTGRU  F2,EDT) 

NOSORT 

*  G.-AR  STIFFMFSS  VALUES  DETERMINED 
RTTH=TH2/AT TH 

M  =  RT  TH 

KAPa(HTTH-M) 

KAP|=1.-KAP 

K34  =  AFGFM( STr  NS.KAP I ) 

K34= 1 . 0F+04-K34 
lF(GU0i)-2)4H,9H,49 
99  K?3  =  AFGPNC  STEMS,  KAP  ) 

K23=K23*1.0F+06 
FL  T  =0. 0 
GD  TO  112 
9R  RM=M 

RNP=NP 

KAMI =RM/RNP 

KAP4=M/(\ti'  +  1  .  /<•  |  . 

1P(KAP1-KAP4)]0"),  106,  106 

105  „23  =  AFGEi'l<  STFNS2.KAP) 

K23-K23*  1 . 0P+06 
1FIK23-1  .  177,  77,71 

77  H23=0. 

GO  TO  79 
7rt  R23=7.92 
79  CONTI  NIIF 
G(1  TO  10  7 

106  K23=AFGFN( STFNS.KAP ) 

K  2  3  =  K2  3* I . OF  +04 

GO  TO  111 

107  CONTINUE 
10H  CONTINUE 

*  INTRODUCTION  OF  IMPACT  TOROUP  INTO  TORSIONAL  DYNAMICS 

KAPFLT =0.620+ H I  mow  1 /  (  2.  *A  t  TH ) 

IF<  KAP-. 620)1 11.104,109 

109  1FIKAP-K APFL1  )  1  10, 111.111 

110  FLT=FLTCAL 
GO  TO  112 

111  FL I =0, 

112  CONTINUE 

*  SAMPLFR 

7  1F1N1-AVG1H.H, 1 1 
ft  CONTINUE 

Y1  =  IMPULS<  SI  IMF  I , (NCR ) 

Y2  = ( 1  IMP -ST  IMP )/ INCK 
Y3  =  ZHIiLI>(  Y  1  ,  Y  2  ) 

Y4-V3+0. 1 
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N«Y4 

FRCIN.NI )=ZHOLO( 

Yt.F  I 

12 

CONTINUE 

1  F ( N-NPTS 111,10. 

10 

10 

N 1  =N  M 1 

STIMF«ST!MF+!NCH 

«NPTS 

11  CONTINUE 
SORT 

TERMINAL 


*  M-T  ANALYSIS  AND  AVERAGER 

on  20  i=i,avg 

00  21  J=1,NPTS 
RSPNI J)=ERC< J.  I  I 

21  CONTINUE 

NR  I  1  )  =  1024 
NO  I  M  =  1 
1  s  i  r,N=- 1 

!E0RM=1 

CALL  E0UR2!  KSPN.NR  .ND1M,  I  SION,  1  f  OR  M  ) 

DO  22  J= 1 , NPTS 

DATA  I  ,I)=0A1  A  1 .1  MR  SEN!  J  l/AVG 

22  CON'INIIF 
20  CONTI  NIJF 

00  2 1  1  =  1 ,  I  OAT  A 
AMPI I ) =CABS( DATA! I  1  ) 

OBAMP( I )=20.»ALO01O( AMP( | l/l.OF-04) 

WI  I)  =  ( I-l, ) / <  |NCR*NPTS) 

* 

*  AMP ( N 1  IS  THF  ABSOLUTE  AMPt  I1U0E  OF  THE  FOURIER 

*  COEFFICIENT. 

* 

23  WRITE!?. 1031WI I l.DHAMPI I) 

103  FORMAT! FR.  1 , F V. ? ) 

00  24  !=1,NPIS 
500  TRA  =  RFAL(FRC(  1,3)1 

TM=(  I-l  )■>  1NCR 
WR  I  TF ( 2 , IOOOITm.TRA 
1000  FORMAT ! 2F 14. 5 1 

24  CONTINUE 

* 

*  I  NCR  IS  THE  SAMPLING  INTERVAL. 

*  INCR  =  1./(2.*F1NT  1  WHERE  PINT  IS  THE  HIGHEST  FRFOUENf, Y 

*  OF  INTEREST  IN  HER  I  Z . 

*  STIMF  IS  THE  TIME  AT  WHICH  THF  ANALYSIS  WILL  START. 

*  CHOOSF  STIMF  SO  THAT  TRANSIENTS  HAVE  OIEO-IHTT. 

*  NPTS  DENOTES  THE  NUMBER  (IE  P0JN1S  IN  THE  FOURIER  TRANSFORM 

*  AVG  IS  THE  PARAMETER  WHICH  DENOTES  THE  NUMBER  OE  AVERAGES 

*  OF  THF  FOURIER  COEFFICIENTS  WHICH  WILL  HF  MADE. 

*  I  DAT  A  DENOTES  THE  NUMBER  OE  FOURIER  COEFFICIENTS  Til  BE  PRINTED. 

* 

PARAM  W0TH=0.02R 
PARAM  NPTS= 1 024 
PARAM  INCR=0. 00005 
PARAM  I  DAT  A  =  500 
PARAM  A VG=4 

PARAM  ST!MF=0.01,ST!ME!=0.01 
PARAM  ISIGN=-1.IEURM=1,N0IM=1 
* 

*  NP  DFNDTFS  THE  NUMBER  OF  TEETH  ON  THF  PINION 

*  CR  DENOTES  THE  CONTACT  RATIO  OF  THE  GEAR  PAIR 

* 

PARAM  FAA=NOOO. 
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PARAM  PJa3. 14159 
PAR AM  NP=4 1 
PARa.”  CRa  1 . 75 

PARAM  J1  a.  010000,  J2  =  .  003903,  J3=.  01  441  ,  J4a.  003?  187,  J5  =  .  0467 
PARAM  B12=6.41,H34al0.HH,R45=l?.04 
PARAM  VPDEL  *3. 95, tfGDFL  =  l. 25 
PARAM  FL  T  AU  =  2  •  5F  -05 
PARAM  K12  =  R»4  3F  +04, K 45  =  5, OF +05 
PARAM  R2  =  1. 92635,  K3  =  2.58<>,R4a  1.92635 
PARAM  FW=0, 2 
* 

*  SSHAFT  DENOTES  INPUT  SHAFT  SPEED  IN  RPM. 

* 

PARAM  SSHAF T  =4500. 

PARAM  KfiaO. 0633, KTal. 100 
PARAM  ROaO. 1 , RL  =0. 1 
PARAM  FLTM=. 05E-03 
PARAM  PS10=-3.9530F-03 
PARAM  PS?D=-1 .7543E-05 
PARAM  PS30a-?. 1959P-05 
PARAM  PS4D=-6.5783F  -04 
PARAM  PS1V=2. 73R7E-01 
PARAM  PS2V=-5.2H71E-01 
PARAM  PS3V=-3. 8600E-01 
PARAM  PS4V=4.4559F-01 
PARAM  TH2l)=1.2575E+01 
PARAM  F 1 =-7 . R500F-02 
PARAM  E2=-1.2263F-05 
* 

* 

*  SFT  DELT  TO  INCR  OR  SUCH  THAT  INCR  IS  A  MULTIPLE 

*  OF  OELT. 

TIMER  OELTaO. 000005 . F  I  NT ! M-0. 32 ,  PROEL =0. 001 , OUTDEL  =  0. 01 

METHOD  RKSFX 

END 

STOP 
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APPENDIX  III 

GEOMETRY  ASSOCIATED  WITH  STIFFNESS  CALCULATIONS 


In  performance  of  the  computation  of  the  stiffness  function  for  the  gears, 
a  set  of  geometrical  relationships  between  the  parameters  used  in  the 
analytical  model  must  be  developed.  To  conveniently  describe  these  rela¬ 
tionships,  a  new  parameter,  t),  is  first  introduced. 

r,  =  distance  from  the  pitch  point  of  the  ^cint  of  contact  between 
mating  teeth  along  line  of  action,  in. 


1  is  positive  for  clockwise  rotation  from  the  pitch  point,  zero  at  the  pitch 
point,  and  negative  for  counterclockwise  rotation  from  the  pitch  point.  To 
begin  the  analysis,  refer  to  Figure  100  and  define  the  following  points: 


A  =  point  of  initiation  of  contact 

P  =  pitch  point 

T  =  point  of  tooth  contact 

B  =  point  of  final  contact 

Oq  =  gear  center 

Op  =  pinion  center 

<t>  =  pressure  angle,  rad 

Rbg  =  radius  of  base  circle  of  gear,  in. 

RpQ  =  radius  of  pitch  circle  of  gear,  in. 

Rtq  =  radius  of  contact  point  of  gear,  in. 

Rgp  =  radius  of  base  circle  of  pinion,  in. 

Rpp  =  radius  of  pitch  circle  of  pinion,  in. 

RTp  =  radius  to  contact  point  of  pinion,  in. 

The  first  parameter  to  be  defined  is  the  radius  of  curvature  to  be  used  in 
the  calculation  of  the  Hertzian  compliance. 

Rrjp  =  !  at| 

-  |ap|  +  tj 

=  ^BP  $  +  TJ  (58) 

Rg  =  |BT| 

-  m  -  b 

=  Rrg  tan  ®  -  T  (59) 

Right  triangle  laws  allow  the  determination  of  the  radius  to  the  contact 
point. 

RTp  =  {(|  AOp  )2  +  (|  AT!  )2|  k 

=  {(rbp)2  +  (Rpp  sin  4  +  (60) 
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100.  Geometry  Associated  With  Gear  Pair. 
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(61) 


rtg  ■  ((|bog|)2+  (  |bt|  )2('a 

55  { ( rbg ) 2  +  (rpg  sin  0  ’  t>)2)’A 

An  approximation  is  used  in  the  expression  for  the  gear  tooth  thickness. 
At  the  contact  point,  by  involutometry  the  following  is  satisfied: 

Rp  cos  $  =  R,p  cos  |T  (62) 

lT  =  involute  angle,  rad 


Using  the  following  relationships,  the  tooth  thickness  at  the  pitch  circle 
is  found: 


Rp  cos  <t> 
R T 


(63) 


where 


where 


I 


T  - 


inv  <t>-inv 


Gg  =  defined  parameter 

tT  =  tooth  thickness  at  the  contact  point,  in. 
tc  =  tooth  thickness  at  the  pitch  circle,  in. 
inv  =  involute  function 


tc  = 


Pc  =  circular  pitch,  in. 


inv  $  =  tan  9  -  <J> 


inv  =  tan 


(64) 

(65) 


(66) 

(67) 

(68) 


The  tooth  thickness  at  any  point  other  than  the  contact  point  is  approxi¬ 
mated  from  the  previous  relationships  in  the  following  manner: 


V  = 


Rp  cos  c 
RB+  Y 


(69) 
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(70) 


?  =  tan 


gs^ 


2X  =  2(Rb  +  y) 


/— 

\2RP 


Gs* 


+  inv  <t>  -  inv  | 


) 


inv  |  =  tan  £  -  | 


(71) 


(72) 


Figure  101  is  used  in  the  description  of  the  angle  (e)  which  specifies  the 
direction  of  the  gear  tooth  centerline  with  respect  to  the  pressure  line. 
From  the  figure, 

C  =  is  the  intersection  of  the  line  of  action  with  the  tooth  center 
CD  =  is  a  perpendicular  to  the  tooth  centerline 


From  Figure  100 , 


LBCD  =  0G 

LACOp  +  LOCD  +  LBCD  =  180° 

LO^D  =  90° 

.-.LACOp  »  90°  -  eG 

Z_BAOp  +  Z_A0pC  +  LAC0p  =  180° 

L  BAOp  =  90° 

.\LA0nC  =  0 
P  G 

LA0j£  =  LAOpT  -  P 


LAOpT  =  tan 


. L  /  Hpp  sin  t>  +  t}\ 
\  RBP  ) 


e  = 


°TP 

2Rrp 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 

(81) 

(82) 


9p  =  tan 


'( 


Rpp  sin  <J>  +  V 


Rr 


UTP 

2R„ 


ABP  /  *^TP 
Figure  102  is  now  used  in  the  determination  of  0p  (for  the  pinion) 

LOqCcB  =  90l  -  0q 


(83) 


(84) 
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From  Figure  100, 


Z_B0gCg  +  £-0gCgB  +  LcGB0p  =  l80° 


LCqBOq  =  90° 


.‘.Z-BOqCq  =  0G 


LBOqT  a  tan 


-l  /  RPG  sin 


From  Figure  101, 


0G  =  tan' 


Z-CG0GT  =  ^tg_ 

2rtg 

-1  /  Rrc  sin  *  -  ’i 


\0jP\  = 


COS  0p 

rBP 
COS  0p 


From  Figure  102, 


ylp  =  tajC  ^  "  RE 


-  ^  cos  ep 


lOtfSl-  ^ 

COS  9r 


^BO 
COS  9q 
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yLG  =  I °GCI  "  RBG 


1 

COS  0q 


(94) 


From  Figure  101, 


1  ct|  =  1  at|  -  |ac| 

=  Rgp  tan  (@p  +  Bp)  -  Rgptan  0p 

(95) 

ii 

I 

=  Rgp  £tan  (0p+  Bp)  -  tan  0pJ 

(96) 

From  Figure  102, 

|CPj  =  |BT|  -  |TCg 

=  Rbg  tan(0G  +  pG)  "  RBGtan  0G  (97) 


hG  *  I  CP| 

=  Rbg  [tan(0G  +  BG)  -  tan  0QJ  (98) 


Two  final  parameters  are  needed  in  the  stiffness  calculation;  i.e.,  the 
values  of  tj  at  the  initial  point  of  contact  and  the  final  point  of  con¬ 
tact  .  From  Figure  100  , 


(Rog  -  W* 


+  Rp-[  sin  4> 


RqP  =  outside  radius  of  pinion,  in. 
Fqq  =  outside  radius  of  gear,  in. 


(99) 


T  ,  =  (Rol  -  rb |)/j  -  Rpp  sin  * 

final 

contact 


(100) 


With  the  geometrical  relationships  developed  in  this  appendix,  the  gear 
mesh  stiffness  functions  may  now  be  expressed  as  a  function  of  the  position 
t)  and  the  physical  dimensions  of  the  gears. 
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APPENDIX  IV 

COMPUTER  PROGRAM:  GEAR  STIFFNESS  CALCULATIONS 


Programs  were  written  in  FORTRAN  to  calculate  the  values  of  the  gear  stiff¬ 
ness  functions.  Separate  programs  were  used  for  gears  in  good  condition 
and  gears  with  faults,  but  the  programming  differences  were  slight.  The 
following  parameters  are  required  in  the  operation  of  either  program: 

ROG  =  outside  gear  radius,  in. 

ROP  =  outside  pinion  radius,  in. 

RBG  =  base  circle  radius  of  gear,  in. 

RBP  =  base  circle  radius  of  pinion,  in. 

RPG  =  pitch  circle  radius  of  gear,  in. 

RPP  =  pitch  circle  radius  of  pinion,  in. 

PI  =  3.14159 

PHI  =  pressure  angle,  rad 

E  =  Young's  modulus,  lb /in? 

G  =  shear  modulus,  lb/in? 

NU  =  Poisson's  ratio 
PC  =  circular  pitch,  in. 

RIM  =  input  speed,  rpm 
FWR  =  input  power,  hp 
W  =  face  width,  in. 

PD  =  diametral  pitch 

PTS  =  number  of  points  of  stiffness  function  to  be  calculated 

The  program  which  calculates  stiffness  for  a  gear  with  a  fault  requires  one 
more  parameter;  i.e., 

WDTH  =  width  of  fault  with  respect  to  the  line  of  action,  in. 

The  calculations  use  a  number  of  subroutines,  as  follows: 


THKNS  =  calculates  tooth  thickness  at  position  off  contact  point 

TKPC  =  calculates  tooth  thickness  at  contact  point 

TRAP  =  performs  a  trapezoidal  numerical  integration  on  the  function: 


(101) 


TRAP2  =  performs  a  trapezoidal  numerical  integration  on  the  function: 


f(y) 


(102) 
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*  A  rr~:;~~-'TT»ij|r'.~:~  :ir  •  imr..:.  ....  ...  i  i . 


TKETA  =  calculates  the  angle  between  the  direction  of  tooth  force 
application  and  perpendicular  to  tooth  centerline. 

The  program  calculates  the  stiffness  functions  based  on  the  equations  pre¬ 
sented  in  the  section  on  Theoretical  Analysis  and  Appendix  III.  Compliance 
values  are  calculated  as  a  function  of  ETA(rj)  per  inch  of  face  width. 

COM(J,l)  =  bending  compliance  of  gear  (j  refers  to  q),  in. /lb 
C0M(J,2)  =  bending  compliance  of  pinion,  in. /lb 
CMHZ(J)  =  Hertzian  compliance,  in. /lb 

The  total  compliance  is  then  calculated. 

COMT(j)  =  tooth  pair  compliance,  in. /lb 

The  program  next  calculates  stiffness  functions  in  the  manner  described 
under  Theoretical  Analysis.  The  results  are  then  formulated  to  produce 
punched  cards  for  use  in  the  torsional  dynamic  computation. 

Listings  of  the  programs  with  parameters  used  in  sample  runs  follow. 


167 


C  PROGRAM  FOR  CALCULATION  OF  STIFFNESS  FUNCT I  ONI  GOOD  GEAR) 
DIMENSION  POSI  100),COMH(  IOO), COM! 100,2),COMT( 100) 
DIMENSION  POSLI  101 ),POS«I  101). STFNS2C 101 1.STFNSI  101) 
DIMENSION  ANGI 101 1 • VAL ( 202  ) 

OIMENS ION  ANG2I 101 ) 

REAL  NU, IOTA,  IOTAG,  IOTAP 

COMMON  PHI, PC 

ROP-2.15 

ROG*2.85 

RBP»3. 8527/2. 

RBG«5. 168/2. 

RPP«2.05 

RPG«5.5/2. 

Pl«3. 14159 
PHI «20.*P I / 180. 

EO0.E+06 
G*ll.5E*06 
NU*. 267 
PC*. 3142 
RPM-4500. 

PMR-15./0.7 

H-0.20 

PTS«50. 

ETA1*-S0RT(  R0G**2-RBG**2 )*RPG*S INI  PHI ) 

£TA2*S0RT I R0P**2-R3P**2 ) -RPP*SI Nl PHI ) 

PO* 10. 0 

PN*PI*COSI PHI )/P0 

ETASS*ETA2-PN 

ETASF*ETA1*PN 

T0*PWR*6. 3E+04/RPM 

FL*TO/RPP 

FO*FL/K 

ETA0-(ETA2-ETA1)/PTS 
C  CALCULATES  BENDING  COMPLIANCE 

1*1 

3  ETA.ETA1 
J*1 

4  CONTINUE 

IF! 1-115,5,6 

5  RP*RPP 
RB*RBP 

RT*SORT I  RB**2  +  ( RP*S INI  PHI  )+ETA)**2) 

GO  TO  7 

6  RP«RPG 
RB*RBG 

RT«S0RTIRB**2+|RP*SINIPHI >-ETA>**2) 

7  CONTINUE 

I0TA*THETA(R8,RP,RT,ETA,  I  > 

VL«RB* 1 1  l./COSI IOTA) 1-1. ) 

88  FORMAT  I  *  ' , 16,  3X ,  1 6,  3X,  E 14.  5,  3X, E 14. 5, 3X , E 14. 5 ) 

CMMM.I  12.*TRAPI  VL.RB.RP 1*1 COSI IOTA)  )**2) /E 

CMSH«< 1.2*TRAP2IYL,RB,RP)*ICOS( IOTA)  )**2)/G 

CMNF*I 1.2*TRAP2(  YL,RB,RP)*I SINI IOTA)  )**2)/E 

COM!  J,  I  )  =  (CMMM+CMNF-fCMSH)/H 

POS I J  )  *ET  A 

ETA=ETA*ETAO 

J*J*1 

0*PTS*2 

IF! J-0)4,29,29 

29  CONTINUE 

IF! 1-1)30,30,31 

30  CONTINUE 
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!■!♦! 

GO  TO  3 

C  CALCULATES  HERTZ  IAN  COMPLIANCE 
31  ETA*ETA1 
J«1 

34  CONTINUE 

RTP«SQRT<RBP**2+<RPP*SlNIPHl )*ETAI**2) 
RTG»SORT(RHG**2+|RPG*SlN<PHl )-ETA)**2> 

RCVG*RBG*T AN(  PHI  )-ETA 
RCVP*R8P*TAN(  PHI l-fCTA 
BETAG*TKPCIRPG,RTGI/(2.*RTG> 

BETAP*TKPC(RPP.KTPI/(  2.*RTP) 

IOTAG*THETA(R8G,RPG,RTG,ETA,2) 

IOTAP*TH£TA(KBP,RPP,RTP,ETA, 1 ) 

ALPHAG*IOTAG*flETAG 

ALPHAP«IOTAP+8ETAP 

HG*RBG*( TAN  I ALPHAG I-TANI inTAGII 

HP*R6P*(TAN( ALPHAPl-TANI  IOTAP) I 

B«SORT<  <4.*FO*RCVG*RCVP|*(2.*l 1 . -NU** 2  ))  / 1 E*P I* I RC VG+RCVP) ) ) 

CMHZ*2. * ( 1. -NU**2 I*  (  ALOGI  4.  *HG*HP/B**2  )-NU/ { l.-NU) ) / ( P |*E ) 

COMHI J)*CMHZ/W 

ETA*ETA+ETAO 

J«J*l 

I F ( J-0 1  34, 39,  39 

39  CONTINUE 
L-PTS+1 

WR I TE ( 6,97 ) 

C  CALCULATES  TOTAL  TOOTH  PAIR  COMPLIANCE 
00  40  J«1,L 

COMTI J).( COMHI JI+COMI J,  U+COMI J,2) ) 

WRITEI6,99)J,P0S< J),  COMTI  J), COMHI  J), COM! J.ll.COMl  J, 2) 

40  CONTINUE 
WRITE(6,98!ETAl,i:T42,TQ 

C  CALCULATES  POSITION  OF  AOJACENT  TOUTH  PAIR  COMPLIANCE  WITH 
C  RESPECT  TO  TOOTH  PAIR  OF  INTEREST 

DO  41  1*1, L 
STFNSI  I)*l. /COMTI  I  1 
STFNS2I I )*STFNS(  I  I 
POSLI  I  l*POS( I l-PN 
POSRI  I  )=POS! I  I *PN 

41  CONTINUE 

42  1*1 
41  J-l 

C  CALCULATES  EFFECT  OF  INITIAL  TWO  TOOTH  PAIR  CONTACT 
44  K*J-M 

IF ( K-L  147,47,51 

47  IFIPOSI  II.LT. POSLI JIIGO  TO  49 
IFCPOSI II.GT.POSLIKJ IGO  TO  50 
OPOS*POSL I K l-POSL  I  J ) 

POS1*POS( I l-POSL  (  J I 
P0S2*P0SL( K l-POSI  I  I 

STFNS2I I )=STFNSI  I  I ♦  (  POS 1* STFNS I  J ) +P0S2* STFNS I K ) I /OPUS 
IFIPOSI I l-ET  ASS )48 ,48 , 5 1 

48  CONTINUE 

49  I *1*1 

GO  TO  43 

50  J=J*1 

GO  TO  44 

51  CONTINUE 

52  1*1 

C  CALCULATES  EFFECT  OF  SECOND  TWO  TOOTH  PAIR  CONTACT 

53  J*l 
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54  K*J  +  1 

IFIK-L (57,57,61 

57  IFIPOSI I I.LT.PUSRI J) »G0  10  59 
IFtP0Sm.GT.PUSK(K))r,0  TO  60 
DPOS*POSRIK)-PUSRI J) 

P0S1*P0S( I )-POSH( J» 

P0S2«P0SR(K)-P0S!  I  ) 

STFNS2I I )=STFNS( I ) ♦ ( POS 1* ST F NS t J ) +P0S2*STFNS I K  )) /OPOS 
IFI POSH  )-ETA2>58,58,61 

58  CONTINUE 

59  1*1+1 

GO  TO  53 

60  J«J+1 

GO  TO  54 

61  CONTINUE 

00  70  1*1,1 

WRITE (6, 99 ( I ,POS< I »,STFNS<  I ) , STFNS2I  I  I, POSH  I I.POSRIII 

70  CONTINUE 

C  PREPARES  STIFFNESS  FUNCTION  FOR  FUNCTION  GENERATOR 
AB»S0RT1R0P**2-R8P**2) 

AAS=AB-<  ETA2-ETA1 ) 

AA  SOA* A  T  AN ( AAS/R8P ) 

ABOA*ATAN( AB/KBP I 
ATRVL*ABOA-AASOA 
00  71  J* 1 , 5 l 

ANGI J ) * AT  AN ( ( A AS-ET  A l  +  POSI  J I l/RBPI-AASOA 

71  CONTINUE 

00  72  K* 1 , 5  1 

ANG2(K)*( ANGI K) -ANGI 11 ) )/( ANGI 39) -ANGI 11 )  ) 

72  CONTINUE 
K*  1 1 
1*1 

100  CONTINUE 

VALI I )=ANG2(K) 

VALI I +1 )=STFNS2(K)/1.0E+06 
IFI  1-58)101, 102,  102 

101  1*1+2 
K*K+ 1 

GO  TO  100 

102  CONTINUE 

PUNCH94,  VAL  (  1), VAL I  2), VAL  I  3),  VAL  1 4),  VAL  I  5), VAL 1 6) 

1*7 

103  CONTINUE 

PUNCH95, VALI I  ) , VAL I  1  +  1 ) , VALI  I +2 ) , VAL I  I +3 ) , VAL I  I  +4) , VAU 1  +  5 1 , 

1VALI 1+6), VALI 1+7) 

1*1+8 

IFI 1-54)103,104,104 

104  CONTINUE 

PUNCH96, VAL I  55), VAL I  561, VALI 57), VAL I  58) 

00  105  1=1,57,2 

WRITE  I  6, 1000  I  VALI l I, VALI  1  + l ) 

105  CONTINUE 

94  FORMAT! 'AFGEN  STFNS* • .F6.4,  • ,  • , F6.4, • , • , F6.4, • , • , F6.4,  •  ,  » , 

1F6.4, ', ',F6.4, ',...') 

95  FORMAT! 7X.F6. 4,  •  ,  ■ , F6.4, • ,  •  ,  F6.4,  •  ,  • ,  F6.4, * , • ,F6. 4, ' ,  •  , 

1F6.4, »,F6.4,  ', *,F6.4,  •  i 

96  FORMAT (7X,E6, 4,  ' , ' , E6. 4,  ' , •  ,  F6 . 4,  «  ,  •  ,  F6. 4 ) 

97  FORMAT! • l • ,  14X,  'ETA • , 1 1 X,  'C0MT',9X, 'COMHZ', l OX,  'COMP • , 1  OX ,  •  COMG' ) 

98  FORMAT! '0* , »ETA1  =  • ,E 14. 5, 3X, 'El A2= • , E 14. 5, 3X,  •  TORQUE* • ,E  1 4 . 5 ) 

99  FORMAT!'  ', 16, 3X, 5E14.5 ) 

1000  FORMAT!  •  S2E14.5) 

STOP 
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ENO 

FUNCTION  THKNSIRB.RP,  Y) 

C  CALCULATES  TOOTH  THICKNESS 
COMMON  PH  I, PC 
GS«RP*COS<PHI|/<RB+Y) 

ZETA*ATAN( SORT ( l.-GS**2 l/GS I 
TC-PC/2. 

INVZET*TAN(ZETA)-ZETA 
INVPHI *T ANI PHI ) -PH  I 

THKNS  =  2.*(R8*Y)*1 TC/(  2. *RPI* INVPHI -INVZET  ) 

RETURN 

END 

FUNCTION  TKPC(RP.RT) 

C  CALCULATES  TOOTH  THICKNESS  AT  CONTACT  POINT 
COMMON  PH  I, PC 
TC=PC/2. 

GT  *RP*COS ( PHI )/RT 

ZET  AT  =  ATAN( SORT ( 1 . -GT** 2  I /GT ) 

INVPHI zT  ANI PHI l-PHI 

INVZET*TAN(ZETAT)-ZETAT 

TKPC*2.*RT*( TC/I2.*RPI*INVPH1-INVZET) 

RETURN 

ENO 

FUNCTION  TRAP! YL,RB,RP) 

C  PERFORMS  INTEGRATION  FOR  COMPLIANCE  DUE  TO  BENDING  MOMENT 

COMMON  PH  I, PC 
F  1*0. 

DELT*YL/200. 

YzOELT 

9  YS*  Y 

Y»Y*OELT 

IF( Y-YL ) 10» 10»  21 

10  X2=< I YL-Y)**2I/( THKNSIRB.RP, Y>)**3 

Xl*( (YL-YSI**2)/( THKNSf  RBfRPfYS) )**3 
FI=Fl*OELT*(Xl+X2)/2. 

TRAP-FI 
GO  TO  9 

21  CONTINUE 
RETURN 

END 

FUNCTION  TRAP2I YL.RB.RP) 

C  PERFORMS  INTEGRATION  FOR  COMPLIANCE  DUE  TO  SHEAR  AND  NORMAL  FORCES 
COMMON  PH  I, PC 
F2*0. 

DELT* YL /200. 

YzOELT 

19  Y$«Y 
Y*Y*OELT 

IFI Y-YL  120,20, 22 

20  X3*1./THKNS(RB,RP,Y) 

X4*1./THKNSIRB,RP,YSI 
F2«F2*0ELT*<  X3  +  X4I/2. 

GO  TO  19 

22  CONTINUE 
TRAP2  *F2 
RETURN 

END 

FUNCTION  THETA(RB,RP,RT,ETA,  I) 

C  CALCULATES  ANGLE  THETA  OF  LUAO  FORCE 
COMMON  PHI, PC 
BETAzTKPCIRP,RT)/(2.*RT) 

IFI  1-1  180,80,62 
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80  THETA«ATAN< (RP*S!N(PHI )+ETA(/RB)-BETA 

GO  TO  64 

62  THETA-ATANI <RP*SINIPHI  I -ETA  ) /RB  I -BETA 

84  CONTINUE 

RETURN 

END 
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C  PROGRAM  FOR  CALCULATION  OF  STIFFNESS  OF  GEAR  WITH  FAULT. 
DIMENSION  POS< 1COI.COMHI 100), COM!  100, 2 ) , COMT ( 1  CO) 
DIMENSION  PUSLI 101  I.PUSRI 101 ) , STFNS2I 101),STFNSI 101) 
DIMENSION  ANGI 101 ), VALI202) 

DIMENSION  ANG2I 101  ) 

REAL  MJ, IOTA, IOTAG, IOTAP 

COMMON  PHI , PC , WOTH 

ROP«2.15 

ROG-2.B5 

R8P»3. 8527/2. 

RBG«5. 168/2. 

RPP-2.05 

RPG«5.5/2. 

PI -3. 14159 
PHI«20.*P1/180. 

E«30.E*06 
G*1 1 , 5F+06 
NU-.287 
PC-.3142 
RPM*4500. 

PWR«15./0.7 
W*0. 20 
PTS«50. 

ETAl«-SQRT(ROG**2-RBG**2)+RPG*SINIPHI) 
ETA2»SORT(ROP**2-RHP«*2l-RPP*SIN(PHl ) 

PD-10. 

PN«P I *COS ( PH  I )/P0 

ETAS$*ETA2-PN 

ETASF«ETA1+PN 

W0TH.0.015 

T0«PWR*6.3E+04/RPM 

FL-TO/RPP 

FO*FL/W 

ETAD*( ETA2-ETA11/PTS 

YPM1N*0.094-WDTH 

YPPLS*0.094*WDTH 

C  CALCULATES  BEN01NG  COMPLIANCE 
1*1 

3  ETA-ETA1 
J«l 

4  CONTINUE 

IF! 1-1)5, 5, 6 

5  RP«RPP 
RB»R8P 

RT«SORT ( RB**2+( RP*SIN( PHI (♦ETA  1**2) 

GO  TO  7 

6  RP*RPG 
RB*RBG 

RT«SQRT (RB**2+(RP*SIN(PHI )-ETA)**2) 

7  CONTINUE 

IOTA=THETA(RB,RP,RT,ETA, I » 

YL»RB*( ( l./COSI  IOTA) 1-1.  ) 

88  FORMAT!'  • ,  1 6 , 3X,  1 6, 3X ,  E 14. 5, 3X ,  F  5, 3X ,  E 14. 5 ) 

CMMM«< 12.*TRAP( YL.RB.RP, I  1*1  COS! IOTA)  )**?)/£ 

CMSH*  (  1. 2* TRAP2! YL.RB.RP, I  I* I  COS! IOTA) >**2)/G 
CMNF«(  1.2*TRAP2( YL,RH,RP, l)*<  SIN! IOTA) )**2)/E 
COM!  J,  I ) =  <  CMMM+CMNF+CMSH ) / W 
POS! J)»ETA 

eta«eta+etad 

J-J+l 

0«PTS*2 

IF! J-D)4,29,29 
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29  CONTINUE 

I F (  1-1130,30,  31 

30  CONTINUE 
!■!♦! 

GO  TO  3 

C  CALCULATES  HERTZIAN  COMPLIANCE 

31  £TA«ETA1 
J«1 

34  CONTINUE 

RTP«SORT|R8P**2*(RPP*SIN(PHI)+ETA)**2l 
RTG«SQRT  <RHG->*2  +  (RPG*SIN!PHl  l-ETA 1**2  I 
RCVG«RBG*TAN( PHI ) -ETA 
RCVP«RBP$T  AN ( PH1I+ETA 

6* SORT) <4.<=F0*RCVG*RCVP)*I 2.*  t l.-NU**2 I ) / ( E*P I  * ( RCVG+kC VP ) >  1 

BETAG*TKPCIRPG,RTGI/(2.*RTGI 

BETAP  =  TKPC(KPP,RTP|/(2.oP.TP) 

IOTAP=THETA< RBP.RPP.RTP.ETA, 1 ) 

IOT AG* THETA! RBG.RPG, R TO, ETA,  2 ) 

ALPHAG=10TAG-*BETAG 

ALPHAP=10TAP+HETAP 

COR*0.0 

YL  P  =  RBP* (  (  l./CUSI  I (JT API  I- 1.  I 
WRITE! 6,2000  I  ETA, RTP.RTG.RC VO,  RCVP 
WRITE (6, 2000  IB, BE  TAG, BE  TAP,  IUTAG, IOTAP 
IE! YLP-YPMIN  1300, 300, 36 

35  IE! YLP-YPPLS  I  36, 300,  300 

36  C0R»2.*W0TH 

corb=ylp-ypmin 

IE(B-CORH)38,30,37 

37  B=8-C0RB 

300  CONTINUE 

HG*RBG*( TAN! ALPHAGI-TAN! IOTAGI l-COR 
HP*RBP*( TAN! ALPHAP  l-TAN!  IUTAP |  ) 

CMHZ»2.*(  l.-NU**2l*(  AL0G<4.eHG*HP/B**?)-NU/t  l.-NUI  >  /  (  H  I  *E  ) 

GO  TO  301 

38  CMHZ=1.0E+64 

301  CONTINUE 
COMH( J)=CMHZ/W 

WRITE(6,2000)ALPHAG,ALPHAP,CnR,  YLG.CORB 
WRITEI6, 2 000)B,HG, HP, C7HZ, ETA 
2000  FORMAT! 5E14.5  I 
E  T  A*ET  A+E  T  AD 
J*J*1 

I E I J-0 134,39,39 

39  CONTINUE 
L=PTS*1 

WR I TE ( 6 , 97 1 

C  CALCULATES  TOTAL  TOOTH  PAIR  COMPLIANCE 
00  40  J  =  1  *  L 

C0MT(JI=(C0MH( JI+COM! J,I)+COM( J,2)| 

WRITE16,99)J,P0S! JI.CUMT! JI.COMHI  JI.COM! J,1),C0M(J,2) 

40  CONTINUE 
WRITE!6,98)ETA1,ETA2,T0 

C  CALCULATES  POSITION  OF  AOJACENT  TOOTH  PAIR  COMPLIANCE  WITH 

C  RESPECT  TO  TOOTH  PAIR  UE  INTEREST 

00  41  l  =  1 ,  L 

IE! COMT  I  I )-l, OE +60)73, 74, 74 

73  STENS!  I  1  =  1. /COMT I  I  I 
GO  TO  75 

74  STENS! I  1=0.0 

75  STENS2!  I )  =  STENSI I  I 
POSL ( I  I =POS <  I  l-PN 
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f 


/ 

f' 


V/  .•» 


POSH ( I 1*P0SI I l+PN 

41  CONTINUE 

42  1  =  1 

43  J*1 

C  CALCULATES  EFFECT  OF  INITIAL  TWO  TOOTH  PAIR  CONTACT 

44  K*  J  + 1 

I F I K-L  14  7, 47 , 5  1 

47  IFIPOSI!!. LT.POSLIJI1GO  TO  49 
IF!  POSI  I  l.GT.POSLt  <1  )G(I  TO  50 
DPOS*POSL(K)-PU$L( J) 

POSl=POS< I l-POSLI J) 

P0S2*P0SL(K l-POSI  I  ) 

STFNS2I I  1  =  STFNS I  I  1  ♦ ( PUS  1 *S TFNS I J 1 +POS2* STFNS I K ) ) /DPOS 
IFIPOSI  I l-ETASS 148.48, 51 

48  CONTINUE 

49  1=1+1 

CO  TO  43 

50  J*J+1 

00  TO  44 

51  CONTINUE 

52  1*1 

C  CALCULATES  EFFECT  OF  SECOND  TWO  TOOTH  PAIR  CONTACT 

53  J*1 

54  K* J+ 1 

IFIK-L 157,57,61 

57  IFIPOSI  I l.LT.POSRIJ) IGU  TU  59 
IFIPOSI  I l.GT.POSRIK 1 IGU  TO  60 
DPOS  =  POSR(K l-POSKI J  1 
P0S1  =  P0S I  I l-PUSKI J  1 
P0S2*PUSR<  K l-POSI I  I 

STFNS2I  I  1  =STFNS  I  I  I  ♦(  PUS  1  *  STFNS  I  J  1 +P0S2*  STFNS  I  K  1 1  /DPI)? 
IFIPOSI  I 1-ETA21  58,58, 61 

58  CONTINUE 

59  1=1+1 

GO  TO  53 

60  J*J+1 

GO  TO  54 

61  CONTINUE 

DO  70  1=1,1 

WRITE! 6,99  ,  I.POSI I  1, STFNS!  I  1,STFNS2I  I  1  ,  POSL I  I  1 ,POSKI  I  ) 

70  CONTINUE 

C  PREPARES  STIFFNESS  FUNCTION  FOR  FUNCTION  GENERATOR 
AB=SQRT(ROP**2-RRP**2) 

AAS*AB-(ETA2-ETA1  1 
AASOA  =  ATAN(  AAS/RBP  ) 

ABOA=ATAN(  AB/RBP) 

00  71  J  =  1 , 5  1 

ANGI J) =ATAN( I  A A S-ET A  1 +  POS I J 1 l/RBPl-AASOA 

71  CONTINUE 

00  72  K  =  1 , 5  l 

ANG2I  K 1=1  ANGI K1 -ANGI  1 1  1  1 /! ANGI 39 J -ANGI  11)1 
72  CONTINUE 
K  =  1 1 
1  =  1 

100  CONTINUE 

VAL I  I  1  =  ANG2 1 K  1 

VALI  I  +  l  1=STFNS2! K  1/1.0E  +  06 

IF!  1-58 1101,102,  102 

101  1=1+2 
K«K+ 1 

GO  TO  100 
102  CONTINUE 
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PUNCH94,VAL<  1  f ,  V  AL  (  2  ) ,  VAL  <  3 ) ,  VAL  I  4 ) ,  VAL  (  5 ) ,  VAL  <  6 1 
1*7 

103  CONTINUE 

PUNCH95.VAL! I» , VAL (  I  +  1 ) . VAL! I *2 ) , VAL ( 1  +  3 ) ,  VAL! I *4 1 , VAL 1 1  +  5 ( , 

1VAL! 1+61, VAl  (1+7) 

1*1+8 

IF( 1-54)103*104, 104 

104  CONTINUE 

PUNCH96,VAL!55),VAl(56),VAL!57),VALI58) 

DO  105  1=1,57,2 

WRITE (6,  1000 ) VAL ( I  ) , VAL ( I  +  1  ) 

105  CONTINUE 

94  FORMAT!  'AEGEN  STFNS2* • , Eft. 4, • , • , T6.4, • ,  •  , F6.4, •  ,  • , F6.4, • , • , 

1F6.4, ', ',Fft.4, ',...') 

95  FORMAT!  7X.F6. 4,  • , • ,F6.4,  • , • , F6.4, • , • , F6.4, • , • , Fft.4, • , • , 

1F6.4, ', *,F6.4, ', '.F6.4, 1 ) 

96  FORMAT! 7X.F6. 4, • , • , F6.4, • , • , F6. 4, • , ' , F6.4 ) 

97  FORMAT!  •  1  • ,  14X ,  • ETA • , 1 IX, 'C0MT',9X, 'COMHZ',  10X, 'COMM* , 10X, •COMG* ) 

98  FORMAT!  ' 0 ' ,  • ET A 1 = ' , E 14. 5 , 3X , <£ T A2= • , E 14. 5 , 3X , • TOROUE=  • , E 14. 5 ) 

99  FORMAT!'  ' , 1 6, 3X , 5E 14. 5 ) 

1000  FORMAT ( 2E 14, 5 ) 

STOP 

END 

FUNCTION  THKNS(R8,RP, Y,  I  ) 

C  CALCULATES  TOOTH  THICKNESS 
COMMON  PH  I , PC , WDTH 
GS*RP*COS (  PH  I  )/(R»+Yl 
ZETA*ATAN<  SORT!  l.-GS**2)/GS) 

TC*PC/2. 

INVZET*TAN(ZETA)-ZETA 
INVPHI *T  AN! PHI  l-PHI 
C0R=0. 0 

YMIN-0.094-W0TH 
YPLS*0.094+WDTH 
IF! 1-11110,110,115 

110  IF  (  Y-YMIN ) 115,115,111 

111  IFIY-YPLS)112,115,115 

112  C0R=2.*W0TH 

115  CONTINUE 

THKNS=2.*!R8+Y)*(TC/( 2. *RP»+ INVPHI- INVZETI-COR 
RETURN 

END 

FUNCTION  TKPC(RP,RT) 

C  CALCULATES  TOOTH  THICKNESS  AT  CONTACT  POINT 
COMMON  PHI, PC 
TC*PC/2. 

GT  *RP*COS ! PH  I )/RT 

ZETAT  =  AT  AN! SORT! 1 . -GT** 2 ) /GT ) 

INVPHI =TAN( PHI  ) —PH  I 

INVZET  =  TAN(  ZETAT  ) -ZETAT 

TKPC=2.*RT*<  TC/( 2 . *RP  I  ♦ I NVPH I - INVZ 6T ) 

RETURN 

END 

FUNCTION  TRAP!  YL,RB,RP, I ) 

C  PERFORMS  INTEGRATION  FUR  COMPLIANCE  DUE  TO  BENDING  MOMENT 
COMMON  PH I  ,  PC , WOTH 
F 1  =0. 

DEL  T  =  YL /200. 

Y=DELT 

9  YS*Y 

Y*Y+DELT 

IF!  Y-YL ) 10, 10,21 
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10 


1 


21 

C 

19 

20 

22 

C 

80 

82 

84 


X2«l ( YL-Y)**2>/<  THKNSIRB.RP, Y, I))**3 
Xl«( (YL-YS)**2)/ I THKNSIRB.RP, YS.  I  )  1**3 
Fl»Fl+DELT*(Xl+X2)/2. 

TRAP*F1 
GO  TO  9 
CONTINUE 
RETURN 

END 

FUNCTION  TRAP2«VL,RB,RP,n 

PERFORMS  INTEGRATION  FOR  COMPLIANCE  DUE  TO  SHEAR  AND  NORMAL  FORCES 
COMMON  PHI.PC.WOTH 
F2=0. 

DEL  T  *YL/200. 

V*DELT 

YS*Y 

Y«Y+DELT 

I F I Y-YL ) 20, 20,  22 
X3*1./THKNS(KB,RP. Yt I  I 
X4»1./THKNS(RB,RP,YS, I  ) 

F2*F2+DEL7*( X3+X41/2. 

GO  TO  19 
CONTINUE 
TR AP2=F2 
RETURN 

END 

FUNCTION  THETA(RB,RP,RT,ETA,  I  ) 

CALCULATES  ANGLE  THETA  UP  LOAD  FORCE 
COMMON  PHI, PC 
BETA=TKPC(RP,RT)/(2.*RT) 

IF  <  1-1)80,80,82 

THETA*ATAN(  <RP*SIN(PHI  ) +ET A  I /RB  ) -BETA 
GO  TO  84 

THETA=ATAN<  <RP*SINIPHI  ) -ETA ) /«B (-BETA 

CONTINUc 

RETURN 

END 
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APPENDIX  V 

TRANSFER  FUNCTION  ANALYSIS  OF  TEST  STAND 


In  order  to  determine  the  structural  dynamic  properties  of  the  test  stand 
system,  a  number  of  transfer  function  tests  were  performed  utilizing  the 
Spectral  Dynanu.es  S.D.  1000  E  Mechanical  Impedance  System  and  an  electro¬ 
magnetic  shaker.  The  driving  motor  was  disconnected  from  the  test  stand 
and  excitation  was  provided  at  the  input  to  the  torquemeter.  The  system 
was  preloaded  by  co^rterweighting  the  generator. 

Figure  103  shows  a  driving  point  impedance  plot  taken  for  torsional  excita¬ 
tion  of  the  test  stand  system.  Evident  on  the  plot  are  a  number  of  reso¬ 
nances  and  antiresonances,  A  test  was  also  performed  in  the  bending  mode 
of  excitation;  a  driving-point  plot  for  this  test  is  shown  in  Figure  104. 
Because  of  the  difficulty  encountered  in  providing  a  purely  torsional 
excitation  force,  with  a  single  linear  stroke,  it  may  be  suggested  that 
when  peaks  (such  as  those  at  305  and  ^75  Hz)  occur  in  both  the  torsional 
and  bending  modes,  they  should  be  ignored,  since  they  are  bending  rather 
than  torsional  resonances.  This  would  thus  suggest  major  torsional  reso¬ 
nances  at  65,  97,  285,  and  660  Hz. 

Transfer  functions  were  also  obtained  for  other  parameters.  Figure  105 
shows  a  plot  of  torquemeter  output  to  driving  point  force  ratio. 

Resonances  are  evident  at  positions  determined  as  torsional  resonances 
by  driving  point  impedance  tests. 

Figure  106  shows  a  plot  of  acceleration  at  the  gear  mesh  to  the  driving 
force.  An  accelerometer  was  mounted  on  the  gear  at  a  meshing  tooth  for 
these  tests.  The  accelerometer  was  positioned  to  monitor  primarily  tor¬ 
sional  rotation.  Figure  107  plots  a  transfer  function  between  this  ac¬ 
celerometer  at  the  mesh  to  an  accelerometer  mounted  outside  the  gearbox 
near  the  input  bearing  position.  This  plot  was  used  to  provide  a  trans¬ 
formation  between  the  rotational  accelerations  generated  by  the  simulation 
and  acceleration  signals  which  could  be  monitored  experimentally. 

On-?  najor  consideration  of  validity  of  results  must  be  realized  in  using 
transfer  function  information.  Since  no  measurement  of  the  coherence 
the  two  signals  was  available,  there  may  be  unreliable  data  points. 
It  ;  •  rather  difficult  to  determine  if  this  is  the  case  unless  more 
sophisticated  testing  tools  are  available. 
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Figure  106.  Acceleration  at  Gear  Mesh  to  Driving  Point  Force 


APFENDDC  VI 

CHARACTERISTIC  ROTATIONAL  FREQUENCIES  OF  TEST  STAND  BEARINGS 


A  number  of  rolling  element  bearings  are  present  in  the  test  stand,  both 
in  the  gearbox  itself  and  in  associated  fixtures.  Calculations  were  per¬ 
formed  to  determine  the  characteristic  rotational  frequencies  associated 
with  the  bearings.  Table  VI  shows  the  dimensions  associated  with  the 
bearings.  The  parameters  for  the  gearbox  input  roller  bearing  are  starred 
since  all  the  bearing  data  was  not  available  and  approximations  had  to  be 
"ude.  The  input  speed  of  4500  rpm  was  used  and  the  frequencies  as  shown 
in  Table  VII  were  calculated.  It  might  be  noted  that  no  significant  amount 
of  signal  was  seen  at  these  frequencies  in  the  experimental  data  that  was 
presented.  Significant  amplitudes  at  these  frequencies  are  seen  when 
faults  are  encountered  on  bearings. 


TABLE  VI.  BEARING  DIMENSIONS 


No.  Balls 

Ball  Diameter 
_ (in.) 

Pitch  Diameter 
_ (in.) _ 

Gearbox 

A.  Output 

14 

0.75 

2.0079 

B.  Idler 

9 

0  50 

2.106 

C .  Input 

i.  Ball 

16 

0.3125 

2.126 

ii .  Roller 

16* 

0.32 

2.224* 

Pillow  Blocks 

A.  Main  Shaft 

9 

0.375 

1.811 

B .  Sliprings 

9 

0.4375 

2.1063 

Sliprings 

12 

0.375 

2.363 

Torquemeter 

11 

0.3125 

_ hm _ 
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TABLE  VII.  BEARING  CHARACTERISTIC  FREQUENCIES  (Hz) 


Cage 

Bass 

Ball 

Bass 

Outer 

Race  Pass 

Inner 
Race  Pass 

Gearbox 

A.  Output 

23  5 

172.8 

328.9 

721.1 

B.  Idler 

21.3 

222.2 

191.9 

311.3 

C .  Input 

i •  Ball 

32.0 

499.2 

511.8 

688.2 

11 .  Roller 

32.2* 

523.2* 

515.7* 

684.3* 

Pilxow  Blocks 

A.  Main  Shaft 

29.7 

346.7 

267.6 

407.4 

B.  Sliprings 

29.7 

345.5 

267.4 

407.6 

Sliprings 

31.5 

460.7 

378.6 

521.4 

Torquemeter 

31.4 

1*45  9 

345.0 

480.0 
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LIST  OF  SYMBOLS 


B53 

B34 


Bt45 

E 


F 

'  0 


G 

hg 

bp 

Ji 
J  r> 

J3 

j4 

J5 


K23 


^34 

^tl3 


Pc 

Pm 

RBG 

rpg 

rtg 


linear  damping  coefficient  of  pinion-idler  gear  mesh,  lb/in.  •  s’1 

linear  damping  coefficient  of  idler-gear  gear  mesh,  lb/in. ‘S-1 

torsional  damping  coefficient  of  input  shafting,  in. -lb /rad’s”1 

torsional  damping  coefficient  of  output  shafting,  in. -lb/rad^s-1 

Young's  modulus,  lb/in.2 

tooth  force  per  unit  facewidth,  lb/in. 

acceleration  due  to  gravity,  in./s2 

shear  modulus,  lb/in.2 

distance  from  contact  point  to  the  gear  centerline  along  the 
angle  of  tooth  loading,  in. 

distance  from  contact  point  to  the  pinion  centerline  along  the 
angle  of  tooth  loading,  in, 

inertia  of  shafting  and  input  gear,  in.-lb«s2 

inertia  of  spur  gear  pinion,  in.-lb»s2 

inertia  of  idler  spur  gear,  in.-lb»s2 

inertia  of  spur  gear,  in.  lb*s2 

inertia  of  generation,  in.-lb*s2 

linear  spring  coefficient  of  pinion-idler  gear  mesh,  lb/in. 

linear  spring  coefficient  of  idler-gear  gear  mesh,  lb/in, 

torsional  spring  coefficient  of  input  shafting,  in. -lb/rad 

torsional  spring  coefficient  of  output  shafting,  in. -lb/rad 

circular  pitch,  in. 

normal  pitch,  in. 

radius  of  base  circle  of  gear,  in. 

radius  of  pitch  circle  of  gear,  in. 

radius  of  base  circle  of  pinion,  in. 
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^*er.rf9F 


0!  ,  02  y  O3  y  04  y  6 5 

'te>  t3>  ^4 

6b 

0 


V 


&H 

tFLT 

6flt 

timp 

0 

l 

1 


radius  of  pitch  circle  of  pinion,  in. 
radius  of  contact  point  of  pinion,  in. 
half -thickness  of  tooth,  in. 

height  above  tooth  root  circle  of  imposed  force,  in. 

angular  rotations  of  inertias,  rad 

relative  angular  rotations  between  inertias,  rad 

tooth  bending  due  to  deflection,  in. 

angle  between  direction  of  tooth  force  application 
and  perpendicular  to  tooth  centerline,  rad 

Fo is son's  ratio 

total  deflection  due  to  Hertzian  force,  in. 

time  gear  spends  in  fault,  s 

fault  width,  rad 

length  of  impact,  s 

pressure  angle,  rad 

involute  angle,  rad 

distance  from  the  pitch  point  to  the  point  of  con¬ 
tact  between  mating  teeth  along  line  of  action,  in. 
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